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FOREWORD 


This  report  represents  one  phase  of  an  analytical  investigation  which,  over¬ 
all,  is  aimed  at  the  development  of  generalized  dynamic  requirements  for  advanced 
vehicle  flight  control  systems. 

The  research  reported  here  was  sponsored  by  the  Flight  Control  Laboratory, 
Aeronautical  Systems  Division,  Air  Force  Systems  Command,  as  part  of  Project 
No.  8219>  Task  No.  821904  It  was  conducted  at  Systems  Technology,  Inc.,  under 
Contract  AT  33(6l6)-596l  .'  The  ASD  project  engineers  have  been  Mr.  George  Xensikls 
end  Mr.  Ronald  0.  Anderson. 

The  authors  wish  to  acknowledge  helpful  discussions  with  Mr.  Malcolm  Abztig 
relating  to  the  axis  systems  proposed,  and  Messrs.  Junichi  Taira  and  Vincent 
Kovace-'rich,  who  performed  several  of  the  developments  in  Chapters  III  and  IV. 
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ABSTRACT 


Conventional  aeronautical  stability  and  control  usage  Is  extended  Into  the 
aerospace  flight  regime  with  particular  reference  to  the  definition  of  operating 
points  and  coordinate  systems.  Implied  and  studied  In  this  extension  are  the 
basic  particle-dynamic  equations  of  motion,  the  magnitudes  of  possible  external 
forces,  and  the  flight  regimes  of  fundamental  Interest.  Simplified  solutions  of 
the  particle-dynamic  equations  of  motion  are  obtained,  and  certain  low-frequency 
modes  common  to  both  particle  and  rigid-body  motion  are  explored. 
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Apofocus,  apogee  (also  used  as  subscript);  constant  of  integration 

Equatorial  radius  of  the  earth 

Polar  radius  of  the  earth 

Areal  velocity  constant,  2(dS/dt) 

Drag  coefficient,  D/(l/2)pU^S 
Lift  coefficient,  nW/(l/2)pU§S 
Center  of  gravity 
Drag  vector;  drag 

Complete  elliptic  integral  of  the  second  kind 
Eccentric  anomaly 

Non-dimensional  altitude  -  equation  76 
Force 

Acceleration  due  to  gravity  at  surface  of  earth  (32.17^0  ft/sec  ) 
Gravity  vector;  local  acceleration  due  to  gravity 
Gravitational  constant  (3-^5  x  10"^  ft^/lb-sec^) 

Perturbed  height  of  vehicle,  perturbed  altitude 

Height  of  vehicle  measured  from  the  surface  of  the  earth,  altitude 
Angular  momentum 

Moments  of  inertia  of  the  earth 


Principal  moments  of  inertia  of  the  vehicle 


Lift  vector,  lift 
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Mass  of  body 

Mass  of  principal  gravitating  body;  molecular  weight 
Semilatus  rectum,  paranieter;  ambient  pressure 
Linear  momentum 
Dynamic  pressure 

Distance  between  centers  of  mass 
Circular  orbit  radius 

2 

Universal  gas  constant;  resultant  aerodynamic  force,  R  =  V  L  +  D 
Semi-major  and  semi-minor  axes  of  ellipse 
Apogee  and  perigee  radii 

Path  distance,  length  of  chord,  displacement 
Laplace  transfonn  variable,  s  =  a  +  jtn 
Perimeter  of  an  ellipse 

Total  area  swept  out  by  the  radius  vector  from  the  gravitating  body 
to  the  orbital  vehicle;  reference  area  for  aerodynamic  coefficients 

Time 

Time  of  perigee  passage;  absolute  temperature 
Thrust  vector,  thrust,  non-dimensional  time 
Orbital  period 

Linear  velocity  along  the  X  (or  X  ')  axis 
Perturbed  velocity  along  the  X  axis 

U/  VgT 

Effective  exhaust  velocity 

Linear  velocity  along  the  Y  (or  Y'  )  axis 

Instantaneous  velocity  of  the  vehicle  observed  in  a  particular  axis 
system,  e.g.,  V  (XYZ) 

Total  velocity  relative  to  inertial  space 
Local  circular  velocity,  a^y/g^/r^ 
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V  Mean  orbital  velocity 

av 

W  Linear  velocity  along  the  Z  (or  Z'  )  axis 

w  Perturbed  velocity  along  the  Z  axis;  argmient  of  perigee 

P  Atmospheric  density  decay  rate  with  height 

7  Longitudinal  flight  path  angle;  also  perturbed  flight  path  angle 

7^  Lateral  flight  path  angle 

r  Ideal  power-off  constant  speed  glide  angle  =  tan  ^ ( -D/l) 

A  Increment  change 

A(s)  Denominator  of  vehicle  transfer  functions;  characteristic  equation 

when  set  equal  to  zero 

Ah  Radial  semi variation.  Ah  =  ae  =  (l/2)(r^  -  rp) 

€  Eccentricity  of  conic  section 

i  Inclination  angle  of  trajectory  plane  to  equatorial  plane 

0,  q)  Perturbed  quantities  of  Euler  angles  0,  and  ♦,  respectively 

%  0,#  Euler  angles  (yaw,  pitch,  and  bank  angles)  relating  XYZ  and  XYZ-axes 
♦  Earth's  gravitational  potential 

X  Latitude  of  the  vehicle 

A  Longitude  of  the  vehicle,  measured  from  the  inertial-^  axis 

Longitude  of  the  ascending  node 
U  Earth ' s  oblateness  parameter 

V  True  anomaly 

p  Mass  density  of  air 

T  Rocket  burning  time 

cUq  Mean  motion 

fi  Angular  velocity  of  vehicle  center  of  mass  about  the  geocenter; 

angular  velocity  of  the  XYZ-axis  system  relative  to  inertial  space 

n  (  )  Angular  velocity  observed  in  a  particular  axis  system,  e.g.,  n(XYZ) 
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T  First  point  of  Aries  (the  ascending  node  of  the  sun  as  seen  from  the 

earth) 

Special  Subscripts: 

a  Aerodynamic 

A  Apogee ,  apof ocus 

b  Burned,  corresponding  to  rocket  fuel  bum  time 

c  Circular  orbit 

d  Disturbances 

f  Final 

i  Initial  value 

n  n  moving  part 

o  Initial  value  or  opierating  point 

P  Perifocus,  perigee,  etc. 

NOTE:  For  summary  of  major  axis  and  reference  systems  see  Table  A-II 
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CHAPTER  I 


INTRODUCTION 


Over  the  last  6o  yeaxs  an  enormous  amount  of  engineering  and  scientific 
effort  has  been  devoted  to  the  study  of  aircraft  dynamics .  With  the  current  ex¬ 
pansion  of  flight  regimes  to  include  extra- atmospheric  phases,  it  becomes  desir¬ 
able  to  extend  familiar  aeronautical  stability  and  control  usage  (Ref.  1  to  5) 
to  include  the  more  general  aerospace  situation.  The  present  report  treats  some 
of  the  initial  topics  which  are  essential  preliminaries  to  this  desirable  direct 
extension. 

Stability  and  control  is  but  o.ie  of  the  three  broad  categories — performance, 
stability  and  control,  and  aeroelasticity — which  together  maice  up  the  field  of 
aircraft  dynamics.  Analytically,  these  classifications  correspond  to  considera¬ 
tion  of  the  aircraft  as  a  particle,  rigid  body,  and  flexible  body.  Overlap  and 
fringe  areas  naturally  exist  between  the  categories;  one  of  the  most  iii5)ortant 
is  that  between  performance  and  stability  and  control.  Because  the  general  rigid- 
body  equations  are  complex  and  nonlinear,  they  are  ordinarily  linearized  about 
"operating  points"  which  define  flight  situations  representative  of  "mean"  con¬ 
ditions.  This  procedural  step  results  in  linear  differential  equations  valid  for 
small  perturbations  about  the  operating  points.  Before  this  can  be  done,  how¬ 
ever,  physically  reasonable  operating  points  must  be  in  hand.  These  are  supplied 
by  the  performance,  or  particle  dynamic  equations.  Consequently,  any  stability 
and  control  problem  implicitly  involves  perfonmnce  considerations  even  when  at¬ 
tention  is  confined  only  to  small  departures  from  operating  point  conditions.  One 
purpose  of  this  report  is  to  define  such  possible  operating  points.  Another,  and 
perhaps  more  important  purpose,  is  to  provide  background  information  on  the  gross 
motions  of  the  vehicle.  With  such  data  available,  the  analyst  is  better  equipped 
to  assess  the  regions  of  validity  for  rigid-  and/or  flexible-body  equations  based 
upon  operating  points  which  have  been  selected  primarily  for  analytical  conven¬ 
ience  . 
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In  the  course  of  developing  the  particle  dynamic  equations  of  aerospace 
vehicles  in  terms  compatible  with  aircraft  stability  and  control,  certain  forces 
and  physical  effects  which  ar^  assumed  negligible  in  aircraft  must  be  taken  into 
account.  Consequently,  a  relatively  large  number  of  new  axis  systems  are  required. 
To  allow  as  much  carryover  between  aircraft  and  aerospace  vehicle  dynamics  as  pos¬ 
sible,  these  axis  systems  must  reduce  to  the  conventional  aerodynamic  case  when 
the  physical  effects  which  caused  their  introduction  are  reduced  to  zero.  An¬ 
other  purpose  of  the  report  has  been  to  define  the  necessary  axes  so  that  such 
a  reduction  is  possible. 

Another  fringe  area  between  performance  and  stability  and  control  involves 
the  motions  and  dynamics  of  the  vehicle  center-of-mass .  Although  aerospace 
vehicles  can  have  significant  coupling  forces  between  translational  and  rotational 
degrees  of  freedom,  this  coupling  is  often  quite  weak  for  the  very  lowest  fre¬ 
quency  modes.  These  modes  are  implicit  in  the  particle  dynamics  formulation. 
Consequently,  the  final  purpose  of  this  report  is  to  elucidate  some  of  the 
characteristic  modes  of  perturbed  motion  occurring  in  the  more  complex  rigid- 
body  situation,  which  are  more  simply  revealed  by  a  study  of  particle  dynamics. 
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CHAPTER  II 


GENERAL  EQJATIONS  OF  MOTION 


As  the  first  step  toward  the  goals  mentioned  in  the  introduction,  the  general 
equations  of  motion  for  an  aerospace  vehicle,  considered  as  a  particle,  shall  he 
developed  in  this  chapter.  Initially  the  equations  are  stated  in  general  vector 
form.  The  various  forces  are  then  resolved  into  an  axis  system  which  can  he  trans¬ 
formed  readily  into  axes  suitable  for  study  of  the  rigid  hody  motions.  The  general 
particle-dynamic  scalar  equations  resulting  from  this  resolution  are  then  special¬ 
ized  into  various  simplified  sets  which  are  considered  in  more  detail  in  subsequent 
chapters . 

By  Newton's  second  law  the  time  rate  of  change  of  linear  momentum  equals  the 
sum  of  all  externally  applied  forces. 


In  rocket-powered  flight  the  thrust  forces,  usually  considered  external,  are  pro¬ 
duced  by  the  expenditure  of  vehicle  mass;  and  the  mass  variation  must  be  consid¬ 
ered  in  determining  the  rate  of  change  of  linear  momentum  (Ref.  6  and  7) •  At  time 
t  let  the  linear  momentum  be 


£;  =  “Y 


where  m  is  the  mass  and  V  the  velocity  of  the  vehicle.  Then,  if  at  time  t  +  dt, 
-dm  is  the  mass  which  has  left  the  vehicle  with  an  effective  exhaust  velocity  , 
relative  to  the  vehicle,  the  linear  momentum  will  be 

£2  =  (m  +  dm)(Y  +  dY)  +  (-dm)(Y  +  y^) 

Here  V  +  dV  and  m  +  dm  are  the  velocity  and  mass  of  the  vehicle  at  time  t  +  dt, 
and  y  +  Vg  is  the  effective  exhaust  velocity  relative  to  inertial  space.  It 
should  be  noted  in  passing  that  the  effective  exhaust  velocity  depends  upon  the 
exit  area,  the  differential  between  exit  and  ambient  pressures,  and  the  exit 
velocity  of  the  mass  leaving  the  vehicle. 


The  Incremental  change  in  momentum  from  time  t  to  time  t  +  dt  is  then 


d£  =  Pg  -  =  (m  +  dm)(V  +  dV)  -  dm(Y  +  -  mY 

=  miy  -  Vgdm  +  dmdV 

Dividing  by  dt,  and  taking  limits,  the  time  rate  of  change  of  momentum  becomes, 


_  dV  dm  „ 

dt  ~  “  dt  ■  ^  dt  - 


dY 


or 


dm 


m  -ST  =  F  +  V-  -rr 
dt  —  — e  dt 


=  F  +  T 


(2) 


where  the  thrust,  T,  is  defined  by  v^(dm/dt). 

Thus,  for  a  particle  where  mass  is  expelled  to  produce  thrust,  the  rate  of 
change  of  momentum  can  be  computed  as  if  the  mass  were  constant,  Euid  the  product 
of  the  change  in  mass  per  unit  time  and  the  relative  velocity  between  the  ex¬ 
hausted  mass  and  the  particle  were  an  external  force.  [Althoijgh  g  =  mV,  the  rate 
of  change  dg/dt  ^  m(dy/dt)  +  y(dm/dt)].  Eq,.  (2)  directly  follows  from  Eq.  (l) 
when  the  thrust  is  developed  by  a  momentum  exchange  other  than  one  directly  in¬ 
volving  the  particle  mass,  as  in  a  propeller  (in  the  present  context  such  a  thrust 
would  constitute  an  external  force).  Thus,  if  the  thrust  force  includes  exhaust 
products,  Eq.  (2)  is  correct  in  general  for  nonrelativistic  (i.e.,  total  mass  is 
constant)  particle  vehicles. 

The  remaining  extemail  forces  on  the  vehicle -particle  are  due  to  gravity,  to 
the  motion  of  the  vehicle  through  the  atmosphere,  and  to  extraneous  disturbances. 
The  gravitationSLl  force  is  represented  by  the  vector  quantity  mg,  \diere  g  is  the 
local  acceleration  due  to  gravity.  The  aerodynamic  force  is  resolved  into  two 
components,  the  lift,  L,  normal  to  the  velocity  vector  and  the  drag,  D,  along  the 
velocity  vector.  An  additionail  force,  F^,  is  used  to  represent  the  sxmimation  of 
all  other,  forces  applied  to  the  vehicle.  These  include  (Ref.  8): 

1 .  Radiation  expulsion  from  the  vehicle 

-Y 

Electromagnetic  radiation  (T-5  x  10  Ib/kw) 

Thermal  radiation  (2.0  x  10“9  Ib/cal/sec) 


k 


2.  Ambient  field  forces 

Gravitation  (from  bodies  other  than  the  principal  source) 
Electromagnetic  fields 

5-  Incident  momentum 

Solar  radiation  (1.8  x  10'"^  Ib/ft^) 

Meteorite  bombardment 
Cosmic  ray  bombardment 

All  of  these  forces  are  either  very  small,  or  have  low  probabilities  of  occur¬ 
rence  (e.g.,  meteorites),  for  aerospace  vehicles. 

Summing  all  the  forces,  the  vector  equation  of  motion  becomes 

dV 

m(t)  ^  =  m(t)£  +  I;  +  D  +  T  +  ^  (3) 

where  the  velocity  vector,  V ,  is  implicitly  referred  to  an  axis  system  which 
is  fixed  relative  to  inertial  space.  Since  the  scalar  equations  to  follow  will 
be  I’eferred  to  an  axis  system,  XYE,  having  an  angular  velocity,  fi  ,  with  respect 
to  inertial  space,  the  acceleration  will  be  given  by 

dV 

^  =  V(XYZ)  +  fl(XYE)  X  Y(XYZ) 

The  time  derivative,  Yonz),  is  the  derivative  of  icm)  obtained  by  treatin,; 
the  XYE  reference  system  as  if  it  were  inertial  space.  The  centripetal  accel¬ 
eration  caused  by  the  curved  path  (relative  to  inertial  space)  of  the  XYE  axis 
system  is  taken  into  account  by  the  vector  cross-product  n  x  ^  .  With  these 
modifications,  the  vector  equation  of  notion  becomes 

Y  +  nxY  =  £  +  ^[L  +  D  +  T  +  P^]  (4) 

where  all  of  the  quantities  are  those  seen  by  an  "observer"  having  a  frame  of 
reference  which  is  rotating  at  an  angular  velocity,  Cl  ,  with  respect  to  Inertial 
space . 

Neither  Eq.  (3)  nor  Eq.  (4)  can  be  solved  in  general.  There  are,  however, 
two  extremely  simple  solutions  to  Eq.  (3)  which  are  applicable  to  high  thrust 
situations;  these  will  be  examined  before  approaching  the  difficulties  of  Eq.  (4). 
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A.  ELEMENTARY  HIGH  THRUST  SOLUTIONS 

The  two  cases  of  interest  involve  situations  where 

(1)  All  external  forces  are  negligible  relative  to  the 
thrust . 

(2)  All  external  forces  except  those  arising  from  a  uniform 
gravitational  field  are  negligible  relative  to  the  thrust. 

Implicit  in  the  statement  of  the  above  assumptions  is  the  further  assumption 

that  the  thrust  is  not  an  external  force  but  is  generated  by  mass  expulsion 

(i.e.,  a  rocket  engine). 

Because  case  (l )  is  a  specialization  of  case  (2)  only  the  latter  need  be 
considered  in  detail.  Then  from  Eq.  (3)  for  g  =  constant,  T  =  Ve  (dm/dt),  and 
all  other  forces  taken  as  zero, 

dV  _  T  dm 

dt  ~  &  ^  ffi(t)  dt 


If  the  effective  exhaust  velocity  is  constant  this  can  be  integrated  to  give 


V(t)  =  Y  +  £t  +  Vg  In 


where  Yqj  Mq,  are  the  Initial  (t  =  O)  and  Y(t),  m(t)  the  present  (t  =  t) 
velocities  and  mass,  respectively;  and  Vg  and  Y  are  oppositely  directed 
(Vg  negative) . 

The  displacement,  can  be  found  if  the  mass  flow  is  constant,  dm/dt  =  -c; 
then,  denoting  the  rocket  burning  lime  by  t,  and  the  final  mass,  m^  -  ct,  by 


I  ds  =  /  Y(t)dt  +  /  Y(t)dt 

-'o  Jj 

=  ^  ^  I  ke  In  ( 


mo  -  ct  /  m^ 

— - - )  dt  +  /  yg  In  —  dt 

®o  I  mo 


^  -  nif 
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The  term  on  the  right  multiplied  by  t  will  always  be  negative,  and  therefore  will 
reduce  the  distance  the  vehicle  travels.  Maximum  distance  will  occur  when  the 
burning  time,  t,  approaches  zero.  Then,  although  the  velocity  eqioation  remains 
unchanged,  the  displacement  approaches 

s  =  ^  +  ]^t  +  i  £t2  -  v^t  in  ^ 

Thus  the  most  efficient  propulsion  method  (to  maximize  the  vehicle  displacement) 
is  one  in  which  the  entire  fuel  mass  is  consumed  instantly  so  that  the  impulse  is 
applied  only  to  the  payload,  and  is  not  expended  in  Imparting  velocity  to  the  mass 
of  unbumed  fuel. 

Eqs.  (5)  and  (6)  may  be  expressed  in  terms  of  an  impulse  I  =  Tt,  or  by  a 

"specific  impulse"  based  upon  the  expelled  weight,  (mo  -  nip)g,  and  defined  as 

Tt 

■^P  (mo  -  mf)g 

by  noting  that 


1 

^  ~  dm/dt 


It  I 

mo-mf  ’’  ”mo-mf  ~  "Sisp 


(7) 


Another  common  term  in  propulsion,  the 
"Characteristic  velocity,"  v^,,  is  defined 
as 

^  ^  -  glep  in  g  (8) 

Physically,  if  there  were  no  gravitational 
field,  the  maximum  velocity  of  the  vehicle 
starting  from  rest  would  be  the  charac¬ 
teristic  velocity.  The  functions  of  mp/mg 
pertinent  to  the  velocity  and  displacement 
equations,  Eqs.  (5)  and  (6),  are  plotted 
in  Figure  1 . 


Figure  1 .  Variation  of  the  Instan¬ 
taneous  Mass  Ratio  and 
its  Logarithm 
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B.  FUBDAMENTAL  COORDINATE  SYSTEMS  AND  TOTAL  ACCELERATION 


The  elementary  case  considered  above  is  indicative  of  typical  simplified 
solutions  to  Eq.  (3)  in  its  present  vector  form.  Eq.  (4)  is  not  even  this  tract¬ 
able.  Therefore,  the  next  major  step  is  to  replace  Eq.  (4)  by  equivalent  scalar 
equations.  Before  this  can  be  done,  a  set  of  reference  axis  systems  must  be 
defined  into  which  the  vector  quantities  may  be  resolved. 

The  selection  of  primary  reference  axes  for  the  study  of  vehicle  particle 
dynamics  generally  reflects  the  types  of  answers  desired  from  the  subsequent 
analytical  treatment  and  the  types  of  forces  to  be  described  mathematically.  In 
the  present  study,  the  main  axis  system  is  chosen  to  have  its  origin  at  the 
particle  Itself,  both  for  analytical  convenience  and  so  it  can  serve  as  one  of  the 
axis  systems  required  in  a  study  of  rigid  body  dynamics.  In  general,  there  are 
several  requirements  and  desires  entering  into  the  coordinate  systems  selected 
for  use: 

a.  The  systems  are  only  Intended  to  be  adequate  for  dynamic 
situations  arising  in  aerospace  vehicle  stability  and 
control.  Thus  only  one  principal  gravitating  body  (the 
earth)  has  important  gravitational  effects  and  the  problem 

.  time  spans  are  relatively  short. 

b.  Both  the  systems  chosen  and  the  notation  used  should 
reduce  to  that  conventional  in  aeronautical  stability  and 
control  as  one  limit. 

c.  The  rotational  sequences  in  going  from  one  axis  system  to 
another  should  be  similar.  In  this  way  only  a  very  few 
basic  transformation  matrices  are  required. 

The  implications  of  these  requirements  and  desires  are  treated  below  and  in  the 
appendix. 

The  physlcail  situation  which  will 
ultimately  be  treated  most  extensively 
in  this  report  is  illustrated  in  Fig¬ 
ure  2.  This  restricted  problem  is 
coplanar,  1  .e . ,  the  velocity  vector, Y , 
is  in  the  XZ  plane,  the  Z  axis  is 
directed  to  the  geocenter,  the  origin 
is  fixed  at  the  particle,  and  the  Y 
axis  is  directed  into  the  plane  of  the 
paper.  The  earth  is  taken  as 
space .  In  normal  aerodynamic 
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and  control  problems  the  XYE  system  shown  becomes  a  reasonable  approximation  to 
inertial  space  since  at  normal  aircraft  speeds  the  rotation  of  XYZrelative  to  the 
earth  (a  ’Taetter"  version  of  inertial  space)  is  negligible.  Instead  of  starting 
with  the  relatively  simple  situation  of  Figure  2,  the  complete  series  of  steps 
required  to  derive  quite  general  equations  which  can  then  be  restricted  to  the 
sln^ilified  case  will  be  taken.  The  ptorpose  of  the  procedure  is  not  to  Induce 
coordinate  psychosis,  but  rather  to  illustrate  some  of  the  assumptions  Implicitly 
involved  in  getting  to  the  simplified  situation  which  will  ultimately  be  studied. 

The  first  problem  is  to  define  inertial  space.  For  aerospace  vehicles  an 
axis  system  with  its  origin  at  the  center  of  the  earth,  and  nonrotating  with  re¬ 
spect  to  the  fixed  stars,  is  suitable.  This  system,  denoted  b.s%  y  ^  ,  is  shown 
in  Figure  3.  The  %  axis  is  directed  along  the  north  polar  axis  of  the  earth,  and 

the  axis  is  directed  to  T,  the  First  Point  of  /ries  (the  sun's  ascending  node 

.  * 

as  seen  from  the  earth) .  The  vehicle  is  oriented  relative  to  inertial  space  via 
a  longitude  angle,  A  (measured  from ),  a  latitude  angle,  and  a  radius  r. 

NORTH  POLE 


ARIES-DIRECTED  GEOCENTRIC  SYSTEM  VEHICLE  LOCATION  RELATIVE  TO 

ARIES-DIRECTED  GEOCENTRIC  SYSTEM 

Figure  3*  Definition  of  Inertial  Axes  for  Aerospace  Vehicles 

An  axis  system  with  its  origin  in  the  body  can  now  be  defined  as  shown 
in  Figure  4.  This  set  of  coordinates,  denoted  asX'Y'Z',  is  the  basic  system 
for  the  formulation  of  the  particle  dynamics  equations  in  their  most  general 
form.  (The  specification  that  X'  be  in  the  direction  of  increasing  latitude 


* 

On  the  celestial  sphere  T  is  the  intersection  of  the  earth's  equator  and  the 
sun's  path  as  the  sun  passes  across  the  equator  from  south  to  north,  i.e., 
the  vernal  equinox. 
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XOC‘) 


PHYSICAL  SITUATION 
xyj  Ames  omecTEO  eeoceNTmc  system 

XYl'  NORTHERLY- AND  GEOCENTRICALLY -OIRECTEO 
SYSTEM  WITH  0RI9IN  IN  VEHICLE 

yC  IN  DIRECTION  OF  INCREASING  LATITUDE  X 

r  IN  DIRECTION  OF  INCREASING  LONGITUDE  A 

«'  DIRECTED  DOWN  GEOCENTRIC  RADIUS 

YW  INTERMEDIATE  AXES 


Figure  4.  Definition  of  Vehicle-Centered  Northerly-  and 
Geocentrically-Directed  Axis  System 


will  he  relieved  later  hy  establishing  an  XYZ  system  with  Z  =  Z'  ,  but  with  X 
oriented  in  the  plane  of  equilibrium  flight. )  The  unit  sphere  development  of 
FigLire  5  shows  the  derivation  of  the  x'y'z'  system  from  the  inertial  reference 
.  This  figure  should  make  the  rather  awkward  choice  of  measuring  longi¬ 
tude  from  the  axis  more  understandable,  since  some  such  artifice  is  neces¬ 
sary  if  the  Z  axis  is  to  be  directed  toward  the  geocenter  (in  common  with 
conventional  stability  and  control  definitions)  and  X'y'z'  is  to  be  readily 
derived  from  %  y.  .  The  directional  cosine  relationships  between  these  two 
systems  are  summarized  in  Figiire  6. 


intermediate  system,  Xy'z" 

DERIVED  FROM  %  ^  ^  BY  ROTATION 
ABOUT  %  THROUGH  LONGITUDE 
ANGLE  A 


VEHICLE-CENTERED  NORTHERLY  -  AND 
GEOCENTRICALLY- DIRECTED  SYSTEM  xVz' 
DERIVED  FROM  XTK  BY  ROTATION  ABOUT  Y 
THROUGH  LATITUDE  ANGLE  X 


Figure  5.  Unit  Sphere  Development  Showing  Derivation  of  Body-Centered 

Radially  Directed  Axes  System  from  Geocentric  Inertial  System 
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The  components  of  the  velocity  vector,  V,  along  X  Y  Z  are  U  ,  V  ,  and  W  . 
These  are  shown  in  Figure  7  together  with  the  definition  of  lateral  and  longi¬ 


tudinal  flight  path  angles. 
The  longitudinal  flight  path 
angle  is  measured  between 
the  velocity  vector,  Y,  and 
its  projection  on  the  X'Y' 
plane,  and  the  lateral  flight 
path  angle  is  measured  from 
the  X'  axis  to  this  same  pro¬ 
jection  of  V.  If  j_' ,  m'  , 
and  n^'  are  defined  as  unit 
vectors  along  X' ,  Y' ,  and 
Z',  respectively,  then  the 
velocity  vector  becomes 


U  (along  X') 


Ycosr  \ 


I  (along  Y') 


Ysin(-y) 


(along  Z  ) 


Figure  7*  Components  of  the  Vehicle  Velocity 
in  the  X'Y'z'  Coordinate  System 


V  =  V(X'y'z')  =  UI'  +  Vm'  +  Wn' 


where 


U  =  Y  cos  y  cos  y-^ 
V  =  y  cos  y  sin  y^ 
W  =  -  V  sin  7 


The  flight  path  angles  are  defined  as 


7  s  sin" 


-W  _  dr/dt 


7i  =  tan-  j 

In  terms  of  components  in  the  X'  Y'  E'  coordinates,  the  angiilar  velocity  of 
the  X'  Y '  Z '  system  relative  to  inertial  space  is 


n(x'Y'z')  =  7[v/-um'] 


and  the  total  acceleration  [the  left  side  of  Eq.  (^4^)3  then  becomes 
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Y  +  ft  X  Y  = 


UJ^'  +V2n'  +  Wn'  +  i 


l'  '  ^1 

1  2!  n  I 

V  -U  0 
U  V  w 


.  (u- if.),'  MV- (w  („) 


The  applied  accelerations  ^the  right  side  of  Eq.  (4)]  are  developed  in  the  follow¬ 
ing  sections. 

C.  RESOU/TION  OF  THE  AERODYNAMEC  FORCES  INTO  X'  Y'  Z'  COORDINATES 

The  aerodynamic  force  acting  on  a  body  is  normally  resolved  into  two  compo¬ 
nents:  the  lift,  L,  normal  to  the  velocity  vector  and  the  drag,  D,  along  the  veloc¬ 
ity  vector  ( see  Figure  8) .  As 
normaiJLy  computed,  lift  and  drag 
depend  upon  the  relative  velocity 
between  the  air  mass  and  the  ve¬ 
hicle,  and  not  the  total  velocity, 

Y .  This  does  not  prevent  resolv¬ 
ing  the  total  aerodynamic  force 
relative  to  the  total  linear  velocity,  as  above,  but  such  a  procedure  can  lead  to 
difficulties  if  care  is  not  taken  when  L  and  D  are  actually  computed. 

Since  the  drag  force  acts  along  the  velocity  vector  it  can  be  resolved  into 
the  X'  Y'  E'  system  using  the  flight  path  angles  already  defined,  i.e., 

D  =  -D[cos  COB  7  j_'  +  sin  7  -  sin  y  n']  (12) 

The  lift,  on  the  other  hand,  is  thus  far  defined  only  as  normal  to  the  velocity 
vector,  so  it  can  conceivably  be  anywhere  in  the  normal  plane .  Consequently  an 
aerodynamic  axis  system,  with  V  and  L  as  two  of  the  axes,  must  be.  related  to  the 
X'  Y'  Z'  system.  This  resolution  can  proceed  as  follows: 


Figure  8.  Directions  of  Lift  and  Drag 

Relative  to  the  Velocity  Vector 


* 

The  entire  problem  can  be  skirted  by  assuming  that  the  air  mass  is  stationary 
relative  to  inertial  space  —  an  assumption  that  throws  away  certain  interest¬ 
ing  second  order  effects  due  to  atmospheric  rotation  caused  by  the  earth's 
rotation. 
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a.  Initially,  let  Y  be  aligned  withX',  with  L  directed  along  -Z' . 

b.  Rotate  Y  about  Z'  through  the  lateral  flight  path  angle, 

(L  is  still  directed  along  -Z'). 

c .  Rotate  Y  about  an  axis  normal  to  both  V  and  Z '  through  the 
longitudinal  flight  path  angle,  y.  At  this  point  Y  is  as 
shown  in  Figure  7. 

d.  Finally,  rotate  L  about  Y,  through  an  aerodynamic  roll  angle,® 
to  its  final  position. 


These  rotations  are  identical  to  those  shown  in  the  Appendix.  (There  a  set  of 
unit  vectors,  _l_,  m,  and  are  related  to  another  set,  i,  ,i,  and  k,  via  rota¬ 
tions  through  'F,  0,  and®.  If  these  angles  are  replaced  by  those  above,  i.e., 

'!'=  y-^,  0=7,  and  ®  =®g,  the  same  form  of  directional  cosine  array  may  be  used.) 
For  the  present  case  this  is: 
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The  lift  force  when  resolved  into  the  X'  Y*  E*  system  is,  accordingly, 

L  =  -l[(cos  y-^  sin  7  cos  ®g  +  sin  7^^  sin 

+  (sin  7-j^  sin  7  cos  ®g  -  cos  7j^  sin  +  (cos  7  cos  ®a)n']  (15) 


D.  GRAVITY  FORCE  EXPRESSED  IN  THE  X'  Y'  E'  SYSTEM 


The  gravitational  acceleration  applied  to  the  vehicle  at  a  distance, 
the  geocenter  woiild  be 


£  = 


from 

(14) 


l4 


if  "the  earth  were  a  homogeneous  sphere*  When  hased  upon  the  arithmetic  mean  of 

the  polar  and  equatoriEil  radii^  the  value  of  g  at  the  earth’s  surface,  the  earth  s 

radius,  a  ,  and  GM  are  about  (Ref.  9) 

^  e 


ge  =  983*05  cm/sec^ 

=  52.252  ft/sec^ 

Eg  =  6.367650  X  108  m 
=  2.089120  X  loT  ft 

GM  =  3.9858  X  10^^  m3/sec2 

=  1  .*10753  X  10*^  ft5/sec2 


These  numbers  are  most  correct  at  latitudes  of  about  *4-5*^  •  Some  idea  of  the  dif¬ 
ferences  due  to  the  actual  earth  shape  and  structure  may  be  obtained  from  the 
table  below,  which  lists  standard  values  for  gravitational  acceleration  of  a  non¬ 
rotating  earth  (Ref.  9) • 


Latitude 

ge 

cm/sec^ 

ft/sec^ 

0° 

981 .43 

52.199 

45° 

983.01 

32.251 

90° 

985.22 

32.258 

To  obtain  a  more  precise  form  of  the  gravitational  force  due  to  the  earth,  a 
better  representation  of  the  earth's  field  is  required.  This  may  be  derived  from 
the  gravitational  potential  function  for  a  homogeneous  ellipsoidal  earth,  with 
moments  of  inertia  I, ^  about  the  polar  axis  and  and  1^^  about  the  principal 
axes  in  the  equatorial  plane.  The  gravitational  potential  can  be  shown  to  be 


*  = 


3  sin^  X) 


3(lyy  ~ 
hVlr^ 


cos^  X  cos  2A  + 
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to  the  first  terms  in  latitude  and  longitude  (Ref.  10).  The  ^  and  ^  axes  used 
■here  are  similar  to  those  defined  in  Figure  5  but  not  quite  the  same;  i.e.,  X  and 
^  are  selected  to  be  principal  axes  of  a  nonrotating  earth,  and  are  not  necessar¬ 
ily  directed  toward  T  .  The  definition  of  longitude  is  also  modified  accordingly. 
Experiments  with  pendiila  indicate  that  the  gravity  vector  varies  only  very  slightly 
with  longitude  so  and  bhe  earth  is  thus  usually  supposed  to  be  approxi¬ 

mately  an  oblate  spheroid.  Under  these  circumstances,  the  second  term  in  the  above 
equation  will  have  the  coefficient 


a 

r 


2 

e 

2 


where  b^  and  a^  are  the  minor  (polar)  and  major  (equatorial)  radii  of  the  earth, 
respectively.  A  close  approximation  to  the  gravitational  potential  of  the  earth, 
including  the  first  order  oblateness  effects,  is  then 


4>  = 


1  +  (i 


(1-3  sin^  X)  + 


(15) 


The  fact  that  the  vehicle  mass  distribution  is  not  actually  that  of  a  point 
introduces  a  small  Increment  in  potential,  ,  given  by 

A«t  =  -  ^1^1  +  3H  (^)^  (1-5  sin2  X)j  di  +  I2  +  I5) 

where  and  I ^  are  the  principal  moments  of  inertia  of  the  vehicle  (Ref.  1 1 ) . 

This  increment  is  entirely  negligible  when  conq)ared  with  Eq.  (15))  so  the  vehicle 
essentially  acts  as  a  point  mass  in  this  instance.  On  the  other  hand,  the  poten¬ 
tial  function  itself  is  not  precise,  since  the  earth  is  probably  pear-shaped, 
nonhomogeneous,  etc. 

To  the  extent  that  Eq.  ( 1 5)  represents  the  gravitational  potential  function 
of  the  vehicle,  the  gravitational  force  can  be  found  by  taking  the  gradient  of 
this  potential. 
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-  grad4> 


Treating  the  coefficients  of  this  expansion  as  enpirical  constants  of  the  earth, 
the  following  sets  of  values  have  been  suggested  (Ref.  12): 

GM  =  3.936277  x  10*^  m^/sec^ 

=  1.i+07732  X  10^°  ft5/sec2 

Eg  =  6.378508  X  10^  m 

=  2.092643  X  lO"^  ft  =  5.963559  X  105  miles 

u  --  5  •  ^6l  X  1 0"^ 

note  that  the  values  of  GM  and  ag  given  here  are  different  from  those  previously 
given  for  a  spherical  earth. 

E.  SCAIAR  FORM  OF  THE  EQUATIOUS  OF  MOTION 

With  the  various  forces  resolved  into  X'  Y'  E'  coordinates,  Eq,.  (4)  can  be 
written  as  the  scalar  set; 
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Along  X'; 


1  i  r  L  /  ^  a\l 

"  5[(ty  -  D  T'l  D  *a  ’'1  ~os  ~)J 
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I D  sin  7  1^1  -  i 
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(IT) 


The  dependent  variables  are  U,  V,  and  W.  Other  quantities  in  the  equations  are 
coimected  to  these  by 


7i  =  tan-1  Y  . 


-W  = 

Q.  V 


7  =  sin"1 


•W 


+  w‘ 


X  =  -ii 


(18) 


L  ®  cj^SC^  p  D  —  (3_SCp 

1  =  2  P^'^rel^  »  ^rel  ”  ^-Yair  mass/inertial  space 


The  atmospheric  density,  p,  is,  among  other  things,  a  function  of  altitude 
and  hence  of  r ;  and  the  thrust  components  may  also  (but  not  necessarily)  be  func¬ 
tions  of  altitude,  speed,  or  dynamic  press\jre.  Thus  all  the  quantities  in  Eqs. 
(17)  are  connected  with  the  three  dependent  variables  in  a  fairly  direct  fashion 
except  the  aerodynamic  roll  angle,  ♦  disturbances,  F;^.  From  the  parti¬ 

cle-dynamics  point  of  view  can  be  an  arbitrary  function  of  time.  Most  of  the 
disturbances  can  also  be  approximated  by  functions  of  time  alone  for  the  periods 


18 


of  interest  here.  Some  of  these  disturbances  have  already  been  mentioned,  but  a 
more  complete  list  is  now  appropriate.  This  list  reflects  assumptions  made  in 
deriving  the  above  equations,  since  is  used  here  as  a  catch-all  to  take  other¬ 
wise  ignored  phenomena  at  least  into  partial  account.  Thus,  to  the  extent  their 
effects  can  be  approximated  for  fairly  short  intervals  by  functions  of  time  only, 
the  con5)0nents  of  might  include  the  following: 

a.  Ambient  Field  Forces 

(1)  Gravitational  effects  of  the  sun,  moon,  and  planets.  These  have 
been  ccmpletely  ignored  in  the  equations  thus  far.  For  close-in 
orbits  and  trajectories,  these  gravitational  attractions  can  be 
approximated  by  periodic  functions  of  time  (with  the  basic  period 
that  of  the  vehicle  about  the  earth) .  When  the  vehicle  is  farther 
out,  or  for  long  time  periods,  this  crude  approach  will  not  suffice. 
In  these  circumstances,  the  whole  formulation  of  the  equations  of 
motion  will  change.  As  a  first  step,  the  "inertial  space"  axes 
would  become  a  heliocentric  ecliptic  system,  the  angular  velocity 
of  X '  Y '  Z '  relative  to  inertial  space  would  be  modified  accord¬ 
ingly,  etc.  Sane  idea  of  when  such  a  transfer  should  be  made  is 
given  by  the  "activity  sphere,"  which  describes  the  space  in  which 
the  planet,  rather  than  the  sun,  should  be  regarded  as  the  center 
body.  Reference  9  cites  an  activity  sphere  of  0.0061 8  astronomical 
units.  This  region  extends  far  beyond  the  "aerospace"  considered 
here.  References  15  and  l4  summarize  the  effects  of  the  sun  and 
moon  upon  close-in  extra-atmospheric  orbits  about  the  earth,  and 
conclude  that  such  effects  are  quite  small.  Over  a  year,  for  exam¬ 
ple,  maximm  changes  in  the  longitude  of  the  ascending  node  are  about 
10  minutes,  and  flucti’ation  in  the  orbital  inclination  is  less  than 
about  4o  seconds.  Consequently,  the  gravitationELL  effects  of  the 
sun,  moon,  and  planets  may  be  ignored  for  almost  einy  conceivable 
aerospace  vehicle  situation. 

(2)  Gravitational  effects  of  the  earth  over  Eind  above  those  represented 
in  the  assumed  gravitational  potential  of  Eq.  (15)~  The  basic  as¬ 
sumptions  leading  to  Eq.  (15)  were  that  the  polar axis  was  a  princi¬ 
pal  SLXis  of  inertia,  that  the  moments  of  inertia  about  the  other  two 
principal  axes  were  equal,  and  that  the  higher  order  oblateness  ef¬ 
fects  were  negligible.  There  is  some  evidence  that  a  difference 
exists  in  the  moments  of  inertia,  but  none  of  the  assiunptions  noted 
are  of  consequence  for  the  time  periods  (several  orbits  as  a  meocl- 
mum)  and  accuracies  of  interest  here . 

(3)  Forces  on  the  vehicle  due  to  electromagnetic  fields.  The  vehicle 
will  normally  have  a  magnetic  moment  and  may  have  a  net  charge.  The 
first  of  these  is  due  to  magnetic  fields  originating  in  intemaO. 
equipment  or,  by  induction,  frcsn  the  earth's  magnetic  field.  The 
net  electrical  charge  could  arise  fran  frictional  effects  within  the 
atmosphere  (normal  bvd.ld-up  of  static  charge),  collection  of  ions, 
etc.  Reference  12,  while  noting  that  a  detailed  analysis  is  very 
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difficult,  also  conjectures  that  the  force  from  these  sources  would 
be  entirely  negligible  relative  to  the  pertiirbation  force  from  the 
oblateness  of  the  earth. 

b .  Incident  Momentum 

(1)  Solar  radiation  forces.  To  be  precise  these  should  be  based  upon 
a  heliocentric  system,  and  the  total  effect  will  be  dependent  upon 
the  actual  configuration. 

(2)  Meteorite  and  cosmic  ray  bombardment. 

c .  Aerodynamic  Disturbances 

Motions  of  the  air  mass  relative  to  inertial  space  enter  into  the 
equations  presently  through  the  lift  and  drag.  Depending  upon  the 
way  in  which  such  air  mass  movements  are  treated,  it  may  be  possible 
to  put  them  into  the  form  of  forcing  functions,  and  hence  into  the 
terms. 

The  above  list  covers  most  of  the  phenomena  which  could  be  taken  into  account,  if 
desired,  as  part  of  the  disturbance  force,  F^.  The  only  additional  terms  that 
might  be  included  are  due  to  radiation  emulsion  from  the  vehicle,  and  these  might 
better  be  considered  part  of  the  thrust  force. 

Besides  the  assumptions  implicit  in  the  above  jtemizations,  Eqs.  (17)  are 
valid  only  \dien  relativistic  and  tidal  friction  effects  are  negligible,  and  when 
the  attitude  motions  of  the  vehicle  are  completely  uncoupled  from  the  translational 
motions.  Reference  12  notes  that  both  tidal  and  relativistic  effects  are  indeed 
negligible  compared  with  other  quantities  included  in  the  equations,  and  that  the 
rotational  and  translational  motions  are  essentially  decoupled  for  extra-atmospheric 
vehicles.  However,  within  the  atmosphere  the  air  forces  can  create  significant 
translational-rotational  couplings . 

F.  SIMPLIFIED  EQUATIONS  OF  MOTION  WITH  NO  OBIATENESS 

The  major  physical  forces  of  interest  for  aerospace  vehicles  are  included  ex¬ 
plicitly  in  Eqs.  (17),  and  the  equations  themselves  are  extremely  complex  and  non¬ 
linear.  The  inclusion  of  the  oblateness  terms,  in  particular,  is  responsible  for 
a  large  amount  of  the  complexity.  The  effects  due  to  oblateness  are  important  for 
long  term  or  extremely  precise  definition  of  close-in  orbits,  but,  as  mi^t  be  ex¬ 
pected  from  the  size  of  n,  these  effects  can  be  safely  ignored  when  only  a  few 
orbital  periods  and  approximate  answers  are  of  interest.  Accordingly  the  first 
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step  in  simplifying  the  equations  for  present  purposes  is  to  discard  oblateness. 
(Oblateness  effects  on  orbital  behavior  will  be  discussed  later  in  this  report.) 

With  the  oblateness  factor  removed,  it  is  no  longer  desirable  to  keep  the 
X'  Y'  Z'  coordinate  system  oriented  as  in  Figure  *4^.  Instead,  it  is  now  possible, 
as  foretold,  to  direct  the  X'  axis  such  that  the  X*  E'  plane  contains  the  veloc¬ 
ity  vector  when  the  vehicle  is  in  equilibrium  flight.  The  vehicle -centered,  geo¬ 
centrically-directed  axis  system  resulting  from  this  step  will  be  denoted  as  XYZ  , 
and  will  henceforth  be  the  reference  system  of  primary  interest.  The  U,  V,  W  and 
7  ,  7 1 ,  notation  for  velocities  and  flight  path  angles  will  be  retained,  al¬ 
though  these  components  are  row  understood  to  lie  along,  or  he  referred  to,  X,  Y, 
and  Z  instead  of  X',  Y' ,  and  Z'.  Inertial  space  may  also  be  redefined  by  an 
^T’  ^T  axis  system.  This  new  version  of  inertial  space  has  its  origin  at 

the  geocenter  and  is  oriented  so  that  the  plane  contains  the  equilibrium 

flight  velocity  vector.  Also  nonrotating  relative  to  the  ^2^^^  Sys¬ 

tem;  and  it  is  not  necessary  at  this  time  to  specify  the  precise  angular  relation¬ 
ships  between  them.  Figure  9  illustrates  the  new  axis  system  structiire,  and  intro¬ 
duces  V  as  the  angular  velocity  of  XYZ  relative  to  inertial  space.  This  is  done 
in  anticipation  of  the  true  anomaly  v  used  in  the  next  chapter. 

With  the  equilibrium  value  of  V  now  contained  in  the  XZ  plane,  the  ratio  V/U 
will  always  be  small  relative  to  unity;  then 

.  V  . 

sin  7^  =  tan  7^^  =  -  ;  cos  7^^  =  1 

and  Eqs.  (17)  sinqilify  to 

^  |-°  7'  [1  '‘’a)]  -  L  sin  ^  +  Tx  +  F^^  j 

VW  1  1  V  V  ) 

T  =  sftT  r°  U  ^  cos4>a  y)  +  Ty  +  F^^^  j  (19) 

=  S  +  j°  sin  7  -  L  cos  7  cos  4>a  +  Tz  +  F^^^  j 
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AXES 

UiriT  VECTORS 

DESCRIPTIOH 

Inertial  system  with  origin  at  geocenter. 
^S^plane  coincident  with  equilibrium 
i^ght  velocity  vector. 

XYZ 

1  m  n 

Rotating  system  with  origin  at  vehicle 
center  of  mass.  Z  directed  toward  geo- 
center.  XZ  plane  coincident  with 
equilibrium  flight  velocity  vector. 

Figure  9*  Modified  Definition  of  Inertial  Space 

and  Radially  Directed  Rotating  Axis  System 
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These  equations  may  be  separated  into  longitudinal  and  lateral  particle 
dynamics  sets  if  the  aerodynamic  roll  angle  is  small  enough,  i.e.,  cos  <I>^  =  1, 
sin  and<I>^V/U  «  sin  7.  The  longitudinal  equations  of  motion  are  then 


•j-  +  5  tan  7)  +  Tx  + 

2 

w+  ^  =  g  +  sin  7  -  L  cos  7  +  Tz  +  F^^j 


(20) 


and  the  lateral  equation,  which  is  still  partially  coupled  to  the  longitudinal 
set  through  U ,  is 

^  =  ir  -  D  cos  7  ^  (1  +  ^  tan  7)  +  Ty  +  F^^]  (21) 

As  noted  earlier,  Eqs.  (20)  could  have  been  simply  derived  by  inspection 
from  Figure  2,  but  all  of  the  implicit  assumptions  involved  would  not  then  have 
been  so  clearly  delineated. 

Although  Eqs.  (20)  and  (21)  ^-fill  be  considered  basic  for  most  of  the  particle 
dynamics  situations  of  this  report,  a  set  that  is  even  simpler  in  come  respects  is 
also  commonly  used.  These  are  equations  in  so-called  relative  wind  axes  which  are 
set  up  by  directing  the  X  axis  along  the  total  velocity  vector  V.  (The  XYE  sys¬ 
tem  would  then  be  true  relative  wind  axes  if  the  atmosphere  were  fixed  relative  to 
inertial  space.)  The  equations  of  motion  for  this  set  of  axes  can  be  derived  from 
Eqs.  (19)  by  suitable  resolution,  remembering  that  the  wind  axis  system  rotates 
with  respect  to  the  radially  directed  system  of  Eqs.  (19)-  A  much  simpler  process 
is  to  write  the  equations  by  insi)ection  of  Figure  9,  using  Eqs.  (19)  as  a  guide 
and  noting  that,  by  definition  of  the  axes,  there  is  never  any  component  of  veloc¬ 
ity  normal  to  Y-  Then, 

Longitudinal : 


-Y7  +  —  cos 


Y  =  -g  sin  7  - 


7  =  g  cos  7 


D 


in(t) 

L 

altT 


+  Tv  +  Fi 


cos  «{>„  +  t4  +  F] 


(22) 
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Lateral: 


vri 


sin*^  +  Ty  +  ^dy 


(23) 


where  the  primed  notation  on  the  T  and  components  indicates  resolution  into 
the  coordinates  based  upon  the  total  velocity,  Y.  Subsidiary  equations,  largely 
definitions,  for  both  Eqs.  (20)  and  (21 )  or  (22)  and  (25)  are 


’'1 

,.2 


dr 

dt 


U  ’ 


y  =  Bin 


2  2 

U  ,* 


.1  rW 


GM 


3  in 


-1  4l/Si  .  t,„-i  d* 


u 


-  W  =  Y  sin  y  =  U  tan  7  =  ^ 


G  = 


1  /  Assumes  air  mass  velocity 

2  ^  2  \ relative  to  inertial  space 

cos  y  ' 


(24) 


Either  set  of  simplified  equations  is  still  so  complex  and  nonlinear  that  no 
general  solution  is  possible.  Therefore,  only  specific  numerical  solutions  can  be 
obtained  if  every  term  in  the  equations  is  of  consequence  to  the  physical  situa¬ 
tion  being  described.  Such  specific  solutions  are  routine  using  automatic  computa¬ 
tion,  but  physical  interpretation  generally  then  requires  calculations  for  a  large 
number  of  cases  involving  suitable  variations  in  fli^t  conditions,  configuration, 
etc.  Accordingly,  it  is  highly  desirable  to  develop  general  solutions  for  further- 
simplified  equations  as  an  adjunct  to  specific  numerical  solutions  of  more  exact 
formulations . 

Further  simplification  can  be  obtained  either  by  neglecting  enough  terms  to 
allow  the  remaining  equations  to  be  solved  readily,  or  by  restricting  the  range  of 
the  variables  to  "small  perturbations"  about  steady-state  "operating  point"  con¬ 
ditions.  When  terms  are  neglected  the  resulting  equations,  which  may  still  be 
nonlinear,  will  have  meaningful  solutions  as  long  as  the  neglected  quantities  are 
physically  negligible .  On  the  other  hand,  keeping  all  the  terms  but  restricting 
their  perturbations  to  "small"  values  always  results  in  a  set  of  linear  differen¬ 
tial  equations.  If  these  are  of  the  constant  coefficient  type  then  solution  is 
simple,  and  the  restrictions  on  the  range  of  the  variables  is  minor  considering 
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the  resulting  general  and  easily  used  information.  When  the  "perturbed"  equations 
are  linear  with  time-varying  coefficients,  no  psirticular  advantage  may  accrue  since 
such  equations  are  often  as  difficult  to  solve  as  the  nonlinear  variety.  Both  c 
techniques  are  useful  devices  for  the  study  of  aerospace  vehicle  peirtlcle  dynamics. 
In  the  two  chapters  immediately  following,  attention  is  confined  to  the  longitudi¬ 
nal  equations  of  motion,  Eqs.  (20),  for  the  two  major  limiting  situations  of  inter¬ 
est  in  aerospace  vehicles.  The  first  corresponds  to  extra-atmospheric  conditions 
where  both  lift  and  drag  are  zero,  and  the  second  to  flight  within  the  atmosphere. 
In  both  instances  power-off  (or  free  flight)  conditions  are  assumed. 
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CHAPTER  III 


SIMPUFIED  EQJATIOMS  FOR  EXTRA -ATMOSPHERIC  FREE  FUGHT 


Extra  atmospheric  free  flight  implies  the  absence  of  lift,  drag,  and  thrust 
forces.  Under  such  conditions  Eqs.  (20)  become 


(25) 


W  + 


r 


(26) 


From  Figure  9  the  angular  velocity  of  XYZ  relative  to  inertial  space,  v,  is 

U 

V  -  7 

^  ••  U  U  .  •  ^  U  . 

and  V  - - -tT  or  U=  rv+— r 

T  fd  r 

Using  these  relationships,  and  noting  that  W  =  ~r  ,  W  =  — r ,  and  g  =  GM/r  , 

Eq.  (25)  reduces  to 

rv  +  2fv  =0  (27) 


r  -  r  = 


(28) 


Eqs.  (27)  and  (28)  are  the  familiar  classical  equations  for  particle  motion  in  an 
inverse  square  field.  The  proi)ertie3  of  this  motion  have  been  well  known  since 
Newton  (e.g..  Ref.  6,  9>  and  15)>  so  only  a  few  salient  features  are  summarized 
below. 


A.  THE  AREAL  VELOCITY  lAW 

The  constant  areal  velocity,  or  angular  momentum,  property  follows 
directly  from  Eq.  (27),  which  is  equivalent  to 

r^v  +  2r?v  =  ^  (r^v)  =  0 

Integrating, 

r^v  =  constant  =  C  (29) 
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Consider  now  an  infinitesimal  element  of  area,  dS,  swept  out  by  the  radius 
vector,  r^,  over  a  path  length  d^,  then 

dS  =  ^  ( r  X  ds) 

Since  d^/dt  is  the  velocity  vector  V,  and  the  linear  momentum  of  the  particle 
is  p  =  mY,  the  areal  velociiy  becomes 


dS  1  ,  ds 
dt  “  2^-^  dt^ 


I  (l  X  Y) 


r_  X  p  is  the  angular  momentum,  H,  in  magnitude  also  equal  to  r  (rv)m  so  the  areal 
velocity  in  terms  of  the  angiilar  momentxnn  per  unit  mass  becomes 


^  ^  H 
dt  m 


C 


(30) 


For  this  reason  the  con¬ 
stant  C,  which  physically 
is  the  angular  momentum 
per  unit  mass,  is  often 
called  the  "areal  veloc¬ 
ity  constant."  Eq.  (30) 
is  a  mathematical  expres¬ 
sion  of  Kepler's  second 
law:  "The  radius  vector 


Figure  10.  Areal  Velocity  Relationship 


between  two  gravitating  masses  sweeps  over  equal  areas  in  equal  times."  This  is 
Illustrated  in  Figure  10. 


B.  INERTIAL  SPACE  VELOCITIES  AND  TRAJECTORIES 


Eqs.  (25)  and  (26)  are  the  components  along  X  and  Z  of  the  vector  equation 
dY/dt  =  -(GM/r^)r^.  In  terms  of  inertial  space  quantities  defined  by;2r^^  coordi¬ 
nates,  the  scalar  equivalent  of  this  vector  equation  (see  Figxire  11 )  is  the  set: 
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GM  , 

=  --^sinv 


Using  the  areal  velocity  law,  =  c/v, 
and  dividing  by  v  =  dv/dt, 


GM  . 

=  -  sin  V 


Separating  the  variables  and  integrating. 


6E0CENTER 


Figure  1 1 .  Geometry  of  Equations  of 
Motion  in  Inertial 
Coordinates 


GM 

%  =  cos  V  +  A 

^  ^  (31) 

/t  =  f  sin  V 

In  Eqs.  (51)  the  inertial  system  has  been  oriented  to  make  one  of  the  constants 

of  integration  zero,  e.g.,  by  directing  the  axes  such  that^^  =  0  when  v  =  0.  When 
this  is  done  the  angle  v  becomes  the  "true  anomaly"  of  celestial  mechanics.  Eqs. 
(31)  can  be  combined  to  obtain  an  expression  containing  both  velocity  con5)onents: 

(52) 

In3^^, coordinates  Eq.  (32)  is  a  circle,  centered  at..2^  =  A,^=  0,  with  radius 

GM/C  as  shown  in  Figure  12.  This  is  the 

hodograph,  or  polar  diagram  of  velocity, 

for  the  motion.  The  toteil  velocity  is 

the  vector  from  the  origin  (pole)  to  the 

circle.  The  "velocity"  of  this  vector 

in  the  hodograph  plane  is  equivalent  to 

the  acceleration  of  the  particle  along 

its  trajectory.  V  andY  .  are  just 
max  min  " 

V  = 

max  C 


Y  =  ^ 

min  C 


Figure  12.  The  Hodograph,  or  Polar 
Diagram  of  Velocity 
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The  trajectory  may  be  found  from  Eqs.  (3I )  without  a  further  Integration  by 
resolving  Into  polar  coordinates, 

0^  =  r  sin  V  ;  HSip  =  f  sin  v  +  fv  cos  v 


=  -  r  cos  V  ; 


r  cos  V  +  r  V  sin  V 


substituting  the  result  into  Eq.  (31 )j  and  eliminating  f  to  obtain 


r  V  =  —  +  A  cos  V 


The  areal  velocity  relationship  is  then  used  to  obtain 


1  GM  ^  A 
r  =  72  ^  C 


which  is  the  expression  for  a  conic  section  in  polar  coordinates  with  the  origin 
at  one  of  the  foci.  This  may  be  more  readily  seen  by  changing  to  the  form 


,  X  AC 

1  ^  '' 


1  +  €  COS  V 


(C^/GM) 


where  e  =  AC/GM  and  p  =  C^/GM  are  the  eccentricity  and  the  parameter  (or  semilatus 
rectum),  respectively.  The  radius  of  the  hodograph,  GM/C,  is  then  the  ratio  of 
the  areal  velocity  constant  to  the  parameter,  C/p,  and  the  center  of  the  hodograph 
is  atiSfcp  =  A  =  Ce/p.  The  maximum  and  minimum  distances  occur  when  v  is  jt  and  zero, 
respectively. 


Ra  =  ,  Rp 


The  corresponding  velocities  are 


V.  .  ^  -  A  . 


Vp  -  f  ‘A 


f  (,  -  O 
f  0 . 0 


which  are  the  minimum  and  maximum  velocities  on  the  holograph.  The  true  anom¬ 
aly  is  thus  seen  to  be  referenced  (i.e.,  v  =  O)  to  conditions  at  perigee;  a 
similar  development  with  the  constanl^  A, changed  in  sign  leads  to  defining  v 
with  respect  to  apogee  conditions.  Both  definitions  have  been  used  in  the 
literature,  but  the  former  seems  to  be  preferred  and  will  be  used  in  this 
report . 

Eq.  (,^)  defines  all  varieties  of  Keplerian  motion.  Some  of  the  fea¬ 
tures  of  the  several  possible  conic  forms  are  summarized  in  Table  I. 

TABLE  I 

VARIETIES  OF  KEPLERIAN  MOTION 


Trajectory 

Eccentricity 

e 

Location  of  Hodograph 
Origin 

IVpI 

|VaI 

Circle 

0 

Center  of  hodograph 

GM 

C 

GM 

C 

Ellipse 

<  1 

Hodograph  includes 
origin 

A 

< 

►u 

A 

Iro 

0  <v^<f 

Parabola 

1 

Hodograph  passes 

2GM 

0 

through  origin 

C 

Hyperbola 

>  1 

Hodograph  excludes 
origin 

>2™ 

For  aerospace  vehicles  the  trajectory  of  primary  interest  is  the 
ellipse;  henceforth,  attention  shall  be  devoted  only  to  this  form.  For 
an  ellipse  the  semimajor  axis,  R^,  is  just  (Ra  +  Rp)/2,  and  the  semiminor  axis 


is 


VTT 


«  Ra^  O'" 


Rn  = 


1  - 

_ p 


(37) 


where 
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Several  of  these  relationships  are  illiistrated  in  Figure  1 5>  which  also  shows 
the  "eccentric  anomaly,"  E,  and  other  astronomical  terms  that  are  used  later, 
including  the  radial  semi  variation,  £:h,  which  is  defined  as 

^  ■  I  (Ra  -  Rp)  =  €Ra  (38) 


By  re-examining  the  several  equations  immediately  above,  it  can  he 
seen  that  any  two  of  the  quantities  e,  R^,  Rp,  R^,  R-jj  and  Ah  may  be  used  to 
describe  the  orbit.  Reference  12  provides  a  very  useful  conversion  table  relat¬ 
ing  these  various  quantities  which,  with  slight  adaptations,  is  reproduced  here 
as  Table  II. 


C.  VELOCITIES  AND  TIME  BEHAVIOR 


Much  of  the  discussion  to  this  point  has  been  concerned  with  the  geometrical 
properties  of  extra-atmospheric  phases  of  aerospace  vehicle  dynamics.  The  veloci¬ 
ties,  both  linear  and  angular,  of  the  vehicle-particle  are  more  important  to  the 
present  report. 

The  total  velocity,  V ,  is  given  by 


2  •2*2 
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TABLE  II 

CONVERSIONS  BETWEEN  GEOMEITRIC  PARAMETERS  OF  EIUIFTIC  ORBIT 
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which,  from  the  hodograph  relationships,  Eq.  (5l)>  can  be  written 


Y 


2 


2GM 


/GM  .  A  \ 

/GM 

A^  \ 

(c2  C  V 

-gm(^ 

■  GM  j 

From  Eq.  (33)  the  first  bracketed  term  is  equal  to  l/r;  and  from  the  definitions 
of  e  and  p  given  in  Eq.  (3^),  the  second  bracketed  term  becomes 

GM  ^  ^  ,  2s 

C2  -  GM  -  c2  ^  ^  ^ 

^  (1  -  £2)  ^  J_ 

P  Ra 

Whereby  Y^  =  GM  (-  -  (39) 


This  is  one  form  of  the  conservation  of  energy  principle,  since  is  proportional 
to  the  kinetic  energy  and  GM(2/r  -  1 /Ra)  to  the  potential  energy  measured  from  a 
constant  reference,  GM/Ra*  At  an  intermediate  point  in  the  above  development  the 
velocity  can  be  obtained  as  a  function  of  the  true  anomaly,  i.e.. 


y2  ^  -h  e  cos  v)  _  1  -  gg] 

L  p  P  J 

,  (^0) 

[Mr 

Y  =  V"p’  ^  ^  ^ 

GM  r  ?  V 

=  [1  +  2€  cos  V  + 

The  quantity  GM/C  which  appears  here  and  in  other  places  throughout  this 
chapter  is,  as  noted  earlier,  the  radius  of  the  hodograph  plot.  It  can  be  delin¬ 
eated  in  physical  terms  directly  from  the  definitions  of  GM  and  C,  but  C  may  it¬ 
self  have  little  direct  physical  significance  for  most  readers.  Accordingly,  a 
few  alternative  expressions  for  this  quantity  'vrtiich  may  lead  to  easier  interpreta¬ 
tion  are  listed  below. 
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GM 

C 


Va  +Vp 
2 


The  last  two  expressions  involve  o^),  the  mean  motion,  and  Vavg>  '•'he  mean  orbital 
velocity,  both  of  which  are  treated  in  the  following  paragraphs  [see  Eqs.  (45) 
and  (47)]*  The  first  expression  follows  directly  from  Eq.  (36),  and  the  second 
is  a  variant  of  the  first.  The  first  is  perhaps  the  most  physically  satisfy¬ 
ing,  and  shows  that  the  quantity  is  the  average  of  the  maximum  and  minimum 
orbital  speeds.  For  a  circular  orbit  with  Va  =  Vp  =  V^,  GM/C  =  Yq  which  follows 
directly  from  either  the  first  three  expressions  above,  or  from  Eq.  (40).  The 
"circular  velocity,"  Vq,  is  discussed  further  later. 

The  flight  path  angle. 


r 


tan 


-_W 

U 


tan 


r 


rv 


tan 


i  IL 

r  dv 


may  be  obtained  by  operating  on  Eq.  (54)  to  give 


tan  7 


€  sin  V 
1  +  e  cos  V 


(^2) 


This  relationship,  converted  into 
the  triangle  of  Figure  l4,  allows 
the  direct  determination  of  the 
horizontal  and  vertical  velocity 
components  U  and  W  as 


Figure  l4.  Flight  Path  Angle  Relationships 
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The  mean  orbital  velocity,  V^yg,  can  now  he  computed  hy  simply  divid¬ 
ing  the  distance  traveled  in  one  revolution  hy  the  period.  The  perimeter  of 
an  ellipse  is 


36 


where  E  is  the  complete  elliptic  integral  of  the  second  kind.  The  mean  orbital 
velocity  is  then 

Yavg  "  T  “  2nRa  ^Ra  “•‘itV^a^  ‘ 

Using  the  approximate  expression  valid  for  ^  1  , 


fGM(2  -  e2)  .  GM 

2Ra  "  C 


f)(l  -  .2) 


V/hcn  the  eccentricity  becomes  zero  the  linear  velocities  all  become 
equal  to  (with  the  exception  of  W,  which  becomes  zero),  i.e.. 


V=  U  =  Yavg  Vc  -  Q 


#1 GM  I—  ISe 

-  y-  .  ViT,  .  .e^- 


where  Y^  is  the  radius  of  the  circular  orbit,  and  gg  is  the  acceleration  due 
to  gravity  at  the  earth's  surface.  The  mean  motion  for  this  same  condition  is 


^  fh  =  X. 
'  If  V  'c  rd  rc  ~  <'c 


The  angular  velocity  of  the  vehicle  on  the  orbit,  and  hence  the  angu¬ 
lar  velocity  of  the  XZ  axes  relative  to  inertial  space,  may  be  obtained  from 
the  areal  velocity  law,  i.e.. 


C  C  ^  .2 

V  =  -5  =  -5  (1  +  e  cos  v) 

r‘^  p<i 


1  -  €‘ 


/ 

/ 1  +  e  cos  v\ 

A 


Separating  the  variables  and  integrating. 


fa^^dt  =  (1  -  €2)5/^r - 

•'  ^  (1  + 


(1  +  G  COS  v)‘ 
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which  becomes,  if  v  =  0  when  t  =  0, 


(jD^t  =  2  tan" 


\1  -  sin  V 
1  +  e  cos  V 


=  2  tan" 


h  -  tan 


=  V  -  tan  7,  when  «  1 

=  V  -  7  (52) 

The  last  step,  Eq-  (52),  follows  from  Eq.  (51)  because  tan  7=7  when  e  is  smaiJ.. 
The  simple  form  of  Eq.  (52)  allows  the  physically  satisfying  picture  of  Figure  15 
to  be  drawn.  Here  the  "mean  anomaly,"  o^t,  is  approximately  given  by  the  angle 

between  the  line  through  perigee  and  a  radial 
line  drawn  normal  to  the  tangential  velocity  vec¬ 
tor  at  the  true  anomaly,  v .  The  relationship 
between  the  mean  and  the  eccentric  ^.nomaly  (see 
Figure  15)  is  given  by  "Kepler's  equation" 

.  (Ref.  15). 


(qjt  =  E  -  £  sin  E 

The  "mean  anomaly,"  ccj^t,  can  also  be 
developed  in  a  series  in  which  cjq-,t  is  a  peri¬ 
odic  function  of  the  true  anomaly.  This  can 
be  accomplished  either  by  expanding  Eq.  (5I ) 
or  by  expanding  the  right  hand  side  of  Eq.  (50) 


Figure  I5.  Approximate  Geo¬ 
metric  Relationship  Between 
Mean  and  True  Anomaly 


before  Integrating.  The  result  (Ref.  15)  is 


cq^t  =  V  -  2e  sin  v  +  — jp 


(' 


sin  2v  -  sin  5v  + 
5 


Another  series,  essentially  the  reverse  of  Eq.  (5^),  which  gives  the  true  anom¬ 
aly  as  a  periodic  function  of  the  mean  anomaly,  can  also  be  derived  (Reference 

15)  • 


a:^jt  +  2e 


5e‘^ 

sin  %t  +  -jj- 


(1  -  sin  2a^t 
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Finally,  the  radius  as  a  function  of  time  (again  with  t  =  0  when  v  =  OJ  can  be 
developed  by  noting  that 


and  from  Eq.  (55) 

1  +  2e  |l  -  -^jcos  fi^t 


(56) 


D.  OKBrr  SPECIFICATIONS 


As  noted  previously,  the  size  and  shape  of  a  particle's  orbit  about 
a  sphere  of  much  larger  mass  is  specified  completely  by  a  knowledge  of  any 
two  of  the  trajectory  parameters: 


Ra 

Rb 

€ 

Ra 

Rp 

zlh 

p 


Semimajor  axis 


Ra>/l 

- 

Semiminor  axis 

f 

Eccentricity 

Rad 

+  e) 

Radius  at  apogee 

Rad 

-  €) 

Radius  at  perigee 

eRa 

1  CM 

< 

II 

Radial  semivariation 

Rad 

-  €2) 

The  parameter,  or  semilatus  rectum 

However,  such  knowledge  is  not  sufficient  to  define  the  spatial  orientation,  as  a 
fiinction  of  time,  of  the  particle  relative  to  coordinates  imbedded  in  the  sphere. 
The  first  additional  requirement  is  to  establish  a  reference  time,  thereby  com¬ 
pleting  the  definition  of  the  trajectory  within  the  orbital  plane.  A  convenient 
reference  time  is  that  corresponding  to  apogee  or  perigee  passage.  Of  these  two, 
perigee  passage  is  the  more  commonly  used.  This  time,  plus  any  two  of  the  above 
parameters,  completely  specifies  the  orbit. 
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Besides  the  location  of  the  particle  on  its  orbit,  the  orientation  of  the 
orbital  plane  relative  to  an  axis  system  in  the  principal  gravitating  body  is  of 
consequence.  In  an  idealized  situation,  a  natural  choice  for  the  reference  axes 
is  a  rectangular  system  fixed  in  inertial  space  with  the  origin  at  the  center  of 
the  gravitating  body.  For  the  case  of  a  spherical  homogeneous  earth  with  no 
oblateness,  the  orientation  of  the  axes  is  a  matter  of  convenience.  The 
axes,  considered  inertial  in  the  last  articles  (and  illustrated  in  Figure  11 ), 
were  chosen  on  this  basis.  When  the  earth's  oblateness  is  to  be  taicen  into  ac¬ 
count,  this  arbitrary  definition  of  inertial  space  is  no  longer  permissible.  In¬ 
stead,  a  basic  reference  frame  like  the  system  shown  in  Figure  3  is  appropri¬ 
ate  .  In  general,  the  orbital  plane  and  line  of  apsides  will  rotate  relative  to 
this  reference  frame.  The  angular  velocities  of  these  rotations  will  ordinarily 
make  only  small  contributions  to  the  total  angular  velocity  of  the  vehicle,  but  it 
is  still  desirable  to  have  some  idea  of  their  magnitude. 

To  present  such  magnitude  information,  and  also  to  develop  the  fact  that  ad¬ 
ditional  quantities  over  and  above  those  already  developed  are  required  to  des¬ 
cribe  the  total  trajectory  situation,  it  is  first  necessary  to  define  angles  and 
axes  which  orient  the  trajectory  plane  relative  to  the  i^^^version  of  inertial 
space.  The  ultimate  axis  set  defining  the  trajectory  plane  is  to  have  the  same 
origin  as  the  .5^!^  inertial  axes  system  —  at  the  focus  of  the  trajectory.  Two 
axes  of  this  set  will  be  chosen  to  lie  in,  and  the  third  to  be  normal  to,  the 
orbital  plane. 

The  trajectory  plane  axes  should  also  be  compatible  with  the  requirements 
and  desires  set  forth  on  page  8.  These  requirements  imply  that  the  axis  normal 
to  the  orbital  plane  be  designated Also,  if  the  coordinate  system  origin 
were  translated  to  the  trajectory,  the axis  should  point  toward  the  geocenter, 
and  the  2^  axis  along  the  trajectory.  The  two  axes  in  the  orbital  plane  should 
also  reflect  the  requirement  that  a  specific  point  of  reference  on  the  trajectory 
is  needed  to  define  a  time  origin.  All  of  these  requirements  are  met  if  the 
axis  is  directed  toward  perigee.  With  the  directions  of^^  and^^  specified  in 
this  way,  the^^  axis  naturally  falls  along  an  appropriate  line  which  is  coinci¬ 
dent  with  the  serailatus  rectum  of  the  trajectory.  The  directional  sense  of  the 
axis  is  not  yet  uniquely  determined,  so  this  sense  is  arbitrarily  made  compati¬ 
ble  with  a  west  to  east  vehicle  trajectory. 
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TRAJECTORY  REFERENCE 
PLANE  GEOCENTRIC  AXES 

Derived  from  %  rotation 

about  f'.  throng  the  co¬ 

inclination;  and  then  about 
through  V,  the  a^unent  of  perigee. 


for/a/  Pton0 

VEHICLE -CENTERED 
GEOCENTRIC  ALLY-OIRECTEO 
AXES  (XY2) 

Derived  fron  by  rotation 

about  ^7  throu^  v,  the  true  anoa- 
aly;  and  then  translation  along 
-B  to  the  vehicle  center-of-mass. 


(a)  (b)  (c) 

Figure  l6.  Orientation  of  the  Trajectory  Plane  Relative  to  Inertial  Space 


This  trajectory  plaxie  axis  system  can  now  be  referred  by  orientation  angles 
to  the  inertial  axis  system  The  evolution  is  illustrated  in  Figure  l6.  As 

the  vehicle  crosses  the  equatorial  plane,  two  orientation  elements  between  the 
equatorial  and  orbital  planes  are  immediately  apparent.  The  first  is  the  point 
at  which  the  vehicle  trajectory  ascends  the  equator,  called  the  ascending  node. 

This  is  located  by  an  angle  Aq ,  the  arc  measured  eastward  along  the  celestial 
equator  from  the  vernal  equinox  to  the  ascending  node.  An  is  referred  to  as  the 
longitude  of  the  ascending  node.  Thus  the  first  intermediate  axis  system  between 
"inertial  space,  and  the  trajectory  plane  axes  is  ascending  node- 

directed  system  of  Figure  l6(a). 

The  second  orientation  element  between  the  two  planes  is  the  orbital  inclin¬ 
ation  angle,  i .  To  take  this  orientation  element  into  account  an  Intermediate 
axis  system  (between  '  and:^^^^j^j  is  formed  by  "rotating about^' 

through  an  angle  n/2  -  i,  the  coinclination.  This  intermediate  system  (not  shown 
in  Figure  l6)  defines  the  orbitail  plane,  i.e.,  but  does  not  provide  a  reference 
point  or  line  (major  axis)  for  the  orbit  itself.  To  do  this,  an  angle,  w,  the 
"argument  of  perigee,"  is  chosen  in  the  orbital  plane  to  measure  the  arc  from  the 
ascending  node  to  the  perigee  point.  The  trajectory  plane  reference  system  is 
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then  established  by  rotating  the  intermediate  axes  about through  w  [see  Fig¬ 
ure  l6(b3- 

With  the  trajectory  reference  plane  geocentric  axes  as  now  developed,  the 
entire  trajectory  is  defined  from  the  viewpoint  of  an  observer  fixed  relative  to 
the  j?i^^-"inertial"  system.  It  should  now  be  apparent  that  six  quantities  are  re¬ 
quired  to  specify  the  total  trajectory  situation,  three  to  define  the  orbital 
plane  and  an  orbit  plane  reference  line,  and  three  to  define  the  actual  vehicle 
time  history  in  the  orbital  plane. 

One  final  axis  system  is  required  to  tie  in  the  present  rather  involved  set 
of  references  with  the  simplified  situation  of  Figures  9>  and  13-  This  is 
the  vehicle-centered,  geocentrically-directed  system  XYZ.  These  axes  are  derived 
from.^^^  by  rotation  about through  the  true  anomaly,  v,  to  form  axes  with 
an  origin  which  is  then  translated  to  the  vehicle  center-of-mass  [Figure  l6(cj]. 
For  normal  aircraft  stability  and  control  purposes,  the  angular  velocity  of  XYZ 
(which  is  just  v  if  the  principal  gravitating  body  is  a  homogeneous  sphere)  is 
ordinarily  small  enough  to  be  neglected.  For  these  circumstances,  the  earth  can 
be  considered  flat,  and  XYE  takes  on  the  character  of  inertial  space. 

When  the  rigid  body  motions  of  the  vehicle  are  considered,  further  axis  sys¬ 
tems  are  required.  Some  of  these  are  illustrated  in  the  appendix,  together  with 
a  summary  of  all  the  systems  discussed  here.  The  appendix  also  contains  figures 
relating  the  Northerly-  and  Geocentrically-Directed  System  (Figure  4)  with  the 
XYZ  axes  via  a  northerly  heading  angle  i'j^. 

For  the  case  of  a  single  central  gravitating  body  which  can  be  considered 
as  either  a  particle  or  a  homogeneous  sphere,  all  of  the  script  axes, 

%'y'^' ,  and  possess  the  qualities  of  inertial  space  in  that  none  of  them 

rotates  relative  to  theiiT^^inertial  reference  frame.  When  spherical  asymmetry 
or  aerodynamic  drag  is  present,  however,  the  orbit  is  no  longer  fixed  relative  to 
inertial  space.  For  example,  oblateness  effects  will  cause  angular  velocities 
Aj^,  di/dt,  w,  and  have  other  Influences  on  the  orbital  situation.  A  cursory  sum¬ 
mary  of  variations  due  to  oblateness  and  drag  is  given  below. 

E.  GENERAL  EFFECTS  OF  OBLATENESS  AND  AERODYNAMIC  DRAG 

ON  THE  IDEAUZED  TRAJECTORY 

As  mentioned  several  times  above,  there  are  a  number  of  forces  which  tend  to 
perturb  the  classical  Keplerlan  orbit  of  a  vehicle.  Of  these,  the  disturbances 
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caused  by  the  attraction  of  the  sun,  moon,  and  planets,  solar  pressure,  the  earth's 
electromagnetic  field,  meteorite  iii5)act,  etc.,  are  all  small  relative  to  those  due 
to  earth's  ohlateness  and  atmospheric  drag.  A  great  deal  of  attention  has  recently 
been  given  to  the  first-order  effects  of  both  these  pertiirbing  influences,  and  the 
qualitative  consequences  are  well  understood.  Quantitatively,  the  effects  of 
oblateness  on  orbital  behavior  are  more  accurately  predictable  than  those  due  to 
air  drag.  This  follows  because  the  magnitude  of  the  oblateness  is  fairly  well 
known,  while  the  detailed  characteristics  of  the  upper  atmosphere  and  the  aero¬ 
dynamic  characteristics  of  vehicles  traveling  therein  are  not.  Also,  to  a  first 
order,  the  oblateness  effects  are  independent  of  vehicle  configuration,  while  aero¬ 
dynamic  effects  are  heavily  dependent  upon  configuration  details . 

In  terms  of  the  five  geometrical  orbital  elements  (s,  w,  l,  and  A^) ,  and 
the  period  (T),  oblateness  and  air  drag  have  the  following  major  first  order  ef¬ 
fects  : 

Oblateness  Effects 

(a)  Secular,  or  aperiodic,  rotation  of  the  line  of  nodes  (causing 
a  A^)  and  of  the  semimajor  axis  relative  to  the  line  of  nodes 
( causing  a  component  of  w) . 

(b)  Periodic  variations  of  the  semimajor  axis  (Rg^)>  eccentricity 
(e),  period  (T),  and  inclination  (i). 

Aerodynamic  Drag  Effects 

(a)  Secular  decreases  in  R^,  e,  and  T,  at  ever  increasing  rates. 

(b)  Periodic  variation  in  the  argument  of  perigee  (w) . 

A  large  number  of  investigators  have  studied  these  effects  for  close  satellites 
(e.g..  Ref.  12  and  l6-22)  .  In  particular.  Ref.  21  presents  a  detailed  treatment 
which  includes  working  equations  for  the  changes  in  orbital  elements  due  to 
oblateness  and  air  drag.  These  changes  are  ejqiressed  there  as  functions  of  the 
true  anomaly,  v,  or  the  eccentric  anomaly,  E,  (both  measiu:ed  from  perigee)  and 
of  the  elements  (Ragj  ^o>  <-o’  "initial"  elliptical  orbit.  This  initial 

elliptical  orbit  may  be  thought  of  as  one  established  about  a  spherical  earth  with 
no  atmosphere.  The  changes  in  orbital  elements  are  then  brought  about  by  "tiirning 
on"  the  oblateness  and  atmospheric  drag  as  if  they  were  thrust  terms.  A  summeu*y 
of  some  of  the  simpler  results  given  in  Ref.  21  is  presented  below  to  give  some 
idea  of  the  magnitudes  and  forms  of  the  various  effects  on  the  orbital  elements. 
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Longitude  of  the  Ascending  Node  (A  )  —  Air  drag  bas  no  effect  on 
the  position  of  the  line  of  nodes.  Oblateness,  however,  rotates  the 
plane  of  the  orbit  about  the  earth's  polar  axis  an  amount,  AAfl,  which, 
for  a  vehicle  in  circular  oibit,  is  given  by 


sin  V  cos  (2wo  +  v)]  cos  Iq 


(57) 


A  similar  result  is  obtained  in  Ref.  12  which  shows  that  a  rotat¬ 
ing  reference  frame  can  be  found  in  which  the  orbital  variations  are 
periodic.  This  reference  orbital  plane  will  have  an  angiilar  rate 


dv 


cos  i 


avg 


Rad  -  e^)J 


cos  i-avg 


(58) 


where  tavg  mean  inclination  of  the  orbital  plane.  The  general 

effect  described  by  Eqs.  (57)  and  (58)  can  be  visualized  as  a  slow  rota¬ 
tion  of  the  line  of  nodes  in  the  earth's  equatorial  plane  as  the  vehicle 
proceeds  on  its  orbital  path.  The  line  of  nodes  will  "regress,"  i.e., 
rotate  clockwise  as  seen  from  the  north  pole,  when  Lq  90°,  and  will 
"progress"  when  Iq  >  90°.  The  existence  of  this  rotating  orbital  plane, 
within  which  the  orbit  is  periodic,  is  one  good  way  of  considering  the 
major  secular,  or  aperiodic,  effects  of  oblateness.  From  Ref.  21 
again,  in  one  revolution  this  rotation  amoimts,  in  general,  to 


and  the  mean  rate  of  rotation  of  the  plane  of  the  orbit  about  the  polar 
axis  is 

.  r  ae 

Afl  =  -10.05  - — ; - ^  COB  i  (deg/day)  (60) 

LRaod  -  eg)J 


This  effect  is  illustrated  in  Figure  17. 


Figure  17.  Regression  of  the  Ascending  Node 
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Inclination  of  the  Orbital  Plane  (i)  —  Drag  does  not  influence  the 
orbital  plane  inclination,  although  oblateness  does  cause  an  oscilla¬ 
tion  of  this  plane  about  the  line  of  nodes.  Uie  amplitvide  of  this 
oscillation  for  a  circular  orbit  is 


Argument  of  Perigee  (w)  —  Both  oblateness  and  air  drag  rotate  the  line  of 
apsides  (major  axis).  The  air  drag,  if  acting  alone,  would  result  in  an  "effective" 
perigee  which,  as  the  vehicle  goes  around  the  earth,  would  oscillate  about  the 
actual  perigee  with  a  period  equal  to  the  orbital  period.  Oblateness  results  in 
a  secular  rotation  of  the  line  of  apsides .  This  rotation  will  be  in  the  same  di¬ 
rection  as  satellite  motion  for  inclination  angles  less  than  about  65-1  /2°,  and 
in  the  opposite  direction  for  higher  values.  The  secular  effect  alone  may  be  con¬ 
sidered  by  examining  the  Aw  which  exists  at  times  when  the  air  drag  effect  is  zero. 
This  occurs  when  the  vehicle  is  at  perigee.  After  one  revolution  the  perigee  has 
moved,  due  to  oblateness,  through  an  angle 

2 

Aw  =  ^ - ^1  (1-5  cos2  Lo)  (62) 


The  mean  rate  for  one  revolution  is  then  obtained  by  dividing  by  the  orbital 
period,  Eq.  (Ml-),  whereby 

.  /  Sp  (l  -  5  cos^  t  ) 

Wavg  =  -5.0  -  (deg/day)  (63) 

A  general  illustration  of  this  secular  effect  is  presented  in  Figure  I8. 


Figure  I8.  Departure  of  Vehicle  from  Conic  and  Rotation  of  Line  of  Apsides 
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Semimajor  Axis  (Ra)  —  Since  the  energy  level  of  an  extra-atmospheric  vehicle 
Is  proportional  to  1 /Ra,  (Eq,.  39) >  any  dissipation  forces,  such  as  air  drag,  will 
be  directly  reflected  In  the  length  of  the  semlmajor  axis.  Is  also  Influenced 
by  oblateness  In  a  periodic  fashion.  This  latter  effect.  Illustrated  In  Figure  l8, 
tends  to  alter  the  conic  form  of  the  trajectory  with  the  vehicle  traveling  sli^tly 
outside  Its  Ideal  elliptical  orbit  when  near  the  equator,  and  slightly  Inside  when 
near  the  apex  of  the  trajectory.  For  a  circular  orbit,  the  change  In  Ra  Is 


ARa 


sin  V  sin  (2wo  +  v)  sln^  lo 


CjjSpRa^ 


m 


(64) 


The  change  due  to  oblateness  is  periodic,  and  vanishes  at  apogee  (v  -  jt)  and  per¬ 
igee  (v  =  O) .  Tlie  seciilar  variation  per  revolution  is  given  by  the  drag  term.  For 
an  assumed  orbital  altitude  of  200  nautical  miles  and  CpS/m  =  6.44,  ARa/Rao  P-^ 
revolution  is  of  the  order  of  2  x  10“^. 

Eccentricity  (e)  —  The  eccentricity  of  -che  orbit  is  continually  decreased  by 
combined  air  drag  and  oblateness,  althou^  the  latter  influence  is  oscillatory  and 
vanishes  at  perigee.  A  primary  effect  of  drag  is  to  "circularize"  the  orbit,  since 
the  apogee  changes  much  more  rapidly  than  the  perigee.  The  oblateness  effect  alone 
is  maximum  at  apogee,  with  a  change  in  "local"  eccentricity  equal  to 


Ae 


+  -y)(l  -  5  sln^  to  Bln^  Wq) 

-  ^  (l  -  Eq)  sin^  to  cos  2Wo 


(65) 


Period  (T)  —  A  decrease  in  period  occurs  due  to  combined  oblateness  and  air 
drag.  For  a  circular  orbit  the  average  change  in  period  per  revolution  is 


AT 

avg 

To 


_  1 

2 


sin'^  tn  cos  2Wr 


^  pSCp 
"2  m 


^o(^  ■  ^e)  (^^) 


Eq.  (66)  indicates  that  there  is  a  difference  in  period  between  a  satellite  on  an 
inclined  orbit  and  one  on  an  equatorial  orbit.-  This  difference,  in  seconds,  con¬ 
sidering  oblateness  alone,  is 


-12.65 


(^^)  '  sln^  to  cos  2wj, 


.1/2 


^ao 


(6?) 


For  more  general  formulas  the  reader  should  see  Ref .  21  . 

It  can  be  appreciated  from  the  foregoing  summary  that,  for  time  spans  of  the 
order  of  a  few  orbital  periods,  all  oblateness  and  drag  effects  are  ordinarily  of 
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negligible  consequence  from  the  standpoint  of  vehicle  stability  and  control  if  any 
reasonable  level  of  automatic  or  manual  control  is  present.  This  statement  is  not 
valid  for  precision  orbit  determination,  many  guidance  situations,  and  very  long 
time,  trimming  only,  flight  control  operations.  In  any  event,  the  above  discus¬ 
sion  and  simplified  formulas  are  helpful  in  determining  whether  or  not  these  ef¬ 
fects  shovild  be  taken  into  account  in  a  given  problem. 
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CHAPTER  IV 

SIMPLIFIED  EQUATIONS  FOR  PCMER-OFF  FLIGHT  WITHIN  THE  ATMOSPHERE 

When  attention  is  confined  to  pov.'er-off  flight  well  within  the  atmosphere 
where  the  air  forces  are  much  larger  than  those  due  to  oblateness  and  external 
disturbances,  a  good  representation  of  the  true  motions  of  the  vehicle -particle 
is  obtained  by  modifying  Eqs.  (20)  and  (24)  to 


To  complete  this  set  requires  that  the  atmospheric  density,  p,  be  connected 
with  the  motion  variables  (Ref-  25)* 


For  a  perfect  gas  the  ambient  pressure  in  terms  of  the  density  is 
given  by 


where  R  is  the  universal  gas  constant,  T  the  ambient  absolute  ten3>er- 
ature,  and  M  the  molecular  weight  of  the  gas.  The  ratio  of  total  dif¬ 
ferential  to  ambient  pressure  is,  then. 


^  ^  ^  ^ 

P  “  p  T  ■  M 


Over  a  large  part  of  the  atmosphere,  especially  those  regions  where 

the  aerodynamic  forces  are  of  major  consequence  in  the  particle  dynamics 

problem, 

^  j  dT  dM 
p  I  ^  T  '  T 
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and 


dp  ^  dp 
P  -  P 

When  the  atmosphere  is  in,  or  approximates,  hydrostatic  equilibriim,  the 
ambient  pressure  variation  with  altitude  is 


dp  =  -pg  dh 

then 

f  - 

and,  to  the  extent  that  p  is  independent  of  h  , 

f 

P  =  Pq®”^^  where  P  =  ^  (69) 

Reference  23  compares  this  exponential  atmosphere  with  the  AEIDC  1 956 
model  atmosphere  when  the  following  values  are  used  for  the  parameters 
listed ; 

Pq  =  0.0027  slug/ft5 

i  =  23,500  ft 

=  50 

Two  figures  appearing  there  are  reproduced  here  as  Figures  19  and  20, 
and  justify,  for  present  piuposes,  the  exponential  model  as  a  first 
approximation  (althoxigh  Pag  will  vary  with  temperature,  etc.,  to  add 
an  additional  ^Q%  or  so  variation  in  Wpag  over  and  above  that  shown 
in  Figure  20) . 

Eq.  (69)  completes  the  set  of  Eqs.  (68),  but  to  obtain  solutions  in  terms  of 
literal  (rather  than  numerical)  values  of  the  parameters  requires  further  simpli¬ 
fication  of  the  equations.  This  can  be  accomplished  either  by  neglecting  terms  or 
by  linearizing  the  equations.  The  first  possibility  is  explored  below,  while 
linearization  is  the  subject  of  the  next  chapter. 

In  the  general  region  of  interest,  the  fractional  change  in  the  radius ,  r  ,  is 
small  relative  to  the  fractional  change  in  forward  velocity,  i.e.. 


or 


(70) 
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DENSITY,^,  tiuf/cii  ft 


Figure  I9.  Comparison  of  Exponential 
Wth  ARDC  1956 
Model  Atmosphere 


Figure  20.  Comparison  of  Mean  Value 
with  AEDC 
Model  Atmosphere 


Above  some  altitude  this  inequality  obviously  will  not  hold,  since  in  extra- 
atmospheric  free  flight  Eq.  (25)  gives 


However,  Reference  25  shows  that  only  a  small  decrease  in  velocity  below  orbi¬ 
tal  conditions  is  required  to  make  the  inequality  of  Eq.  (70)  a  reasonably 
good  one.  For  example,  in  a  typical  re-entry,  dr/r  =  0.1 (du/u)  at  a  point  where 
drag  has  slowed  the  vehicle  by  only  a  percent  or  so  of  its  orbital  velocity. 
Situations  for  which  the  assumption  is  reasonably  valid  are  illustrated  in 
Figure  21  . 

Another  simplification  to  the  equations  stems  from  the  fact  that  kinetic 
energy  changes  are  much  larger  than  potential  energy  changes  for  the  general 
flight  regime  of  interest  here.  The  product  (L/D)  tan  7  is,  accordingly, 
usually  small  relative  to  unity.  This  is  always  true  in  pure  ballistic  (zero 
lift)  flight,  but,  with  lift  present,  the  magnitude  of  the  acceleration,  U, 
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must  be  fairly  large.  In  normal  subsonic  glides,  for  instance,  ij  is  small  and 
(L/d)  tan  7  =  1.  Use  of  the  assunqotion 


Reference  23  treats  these  equations  in  considerable  detail,  and  presents  many 
numerical  solutions  in  terms  of  a  dimensionless  dependent  variable, 

Z  =  pQCj)SU/2m  -^r".  To  obtain  analytical  solutions  in  literal  terms,  however, 
still  further  simplifications  must  be  made;  and  three  relatively  simple  limiting 
cases  are  of  special  interest.  These  are  shown  in  Table  III  with  the  assumptions 
required  (over  and  above  the  two  already  given)  to  reduce  the  equations  to  solv¬ 
able  circumstances. 

As  indicated,  the  simplified  equations  can  be  justified  using  various  as¬ 
sumptions,  some  of  which  are  not  necessarily  as  restrictive  as  others.  The  bas¬ 
ic  assumption,  which  is  also  least  restrictive,  is  given  first.  It  indicates 
the  minimum  condition  under  which  Eqs .  (72)  may  be  truncated  to  those  shown  in 
the  table.  Note  in  particular  that  when  the  third  set  of  assumptions  applies, 
the  "effective"  ballistic  and  the  glide  cases  become  identical.  The  equations 
for  the  sin^jllfied  cases  are  considered  in  more  detail  below.  (See  also  Ref. 

23-27) 
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TABLE  III 

SUMMARY  OF  READILY  SOLVABLE  PARTICLE  DYNAMICS  CASES  WITH  AERODYNAMIC  FORCES 


SITUATION 

ASSUMPTIONS 

EQUATIONS  OF  MOTION 

Effective 

■ballistic 

Gravity,  centripetal,  and  lift  forces 
negligible 

So:  -  g  +  ^  cos  r|  <|w|,  ^  sin  y 

or  g,  ^  cos  r  <C  |w|,  ^  sin  7 

.  L 

or  —  =  g  -  -  cos  7 

•  D 

U  =  -  —  cos  7 
m 

W  =  —  sin  7 
m 

Glide 

Vertical  acceleration  and  vertical 
component  of  drag  negligible 

So:  jw-  ^  sin  7]  <  g>  ^  cos  7 

or  |w|,  J  sin  7  <  g,  J  cos  7 

or  W  =  —  sin  7 

m 

M  D 

U  cos  7 

m 

U2  L 

—  =  g  -  -  cos  7 

Skip 

Gravity  and  centripetal  forces  negli¬ 
gible 

So:  |g  -  ^  •<  Iw],  —  cos  7,  —  sin  7 

1  *  *  131  in 

or  g,  ^  «  |w|,  ^  cos  7,  5  sin  7 

U2  . 

or  —  =  g 

I'l  D 

U  - - cos  7 

m 

W  =  -  —  cos  7  +  —  sin  7 
m  m  ' 
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pia 


Eliminating  D/m, 


But,  by  definition. 


-  tan  y 


W 

U 


dW 
■  dU 


tan  y 


W 

U 


(7^) 


and  for  both  expressions  to  be  valid,  W/U  must  be  constant.  The  simplified 
equations  imply,  therefore,  that  the  flight  path  angle  is  constant,  recognized 
hereafter  by  adding  a  subscript  zero  to  y. 

Eliminating  time  by  dividing  the  first  of  Eqs.  (73)  by  Eq.  (7^) 
and  separating  the  variables. 


u 


Defining  a  nondimensional  altitude,  f,  (Ref.  27) 

Po 


f  s  3h  -  In 


P  sin 


,  df  =  pdh 


_  =  e  df 


Eq*  (75)  becomes 
and,  after  integration,  with  the  initial  velocity  taken  as  Uq, 


SL  _ 

Uo 


Since  dU/dt  =  (dU/df) (df/dt),  the  acceleration,  U,  will  be 


(75) 


(76) 


(77) 


U  =  (Ue"^)(ph)  =  U^p  tan  7oe”^ 

=  PU^  tan  7o  [e'^^ 


(78) 
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The  density  variation  over  the  trajectory  can  be  found  from  the  relationship 


2a  = 

u2 


.2P  (  “  )  e-^ 

^  cos 


Note,  from  Eqs.  (80)  and  (8l),  that  must  be  negative  if  q  and  p  are  to  be 

positive,  as  they  are  in  a  real  physical  case.  Plots  of  normalized  acceleration, 

U/U  ,  and  normalized  velocity,  U/U  ,  are  presented  in  Figure  22  (adapted  from 
niajC  nisix 

Ref.  27). 

2.  Gliding  Flight 

Under  the  restrictions  leading  to  gliding  flight,  the  simplified  eqmtions 
of  motion  are 


D 

U  =  -  -  COB  7 
U2  L 

—  =  g  -  -  cos  7 


(82) 
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Figure  22.  Normalized  Variables  Versus  Altitude,  Ballistic  Flight 


Eliminating  cos  7, 


TI.  Then, 


dt 


dU 


>/^  ^  i  >/i  ^  I 


(85) 
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U  u  2  g  2  r 


From  Newton's  law  of  gravitation, 

GM  dg 

«  =  72  '  T 


dU  _  ^  1  df 
U  ■  u  ■  2  r 


dr 

/  dU 

T 

'  TT 

dU  .  dn  ,,, 

-  =  -  or  dU 


=  -^du  =  \/^  dTI 

u  » 


Recognizing  that  in  the  above  result  is,  by  the  nature  of  the  differenti¬ 

ation  process,  the  initial  rather  than  the  local  value 

dU  4  dtr  =  Y^dn 


(83)  then  becomes 


U  i  YcH 


1  -  u^ 


life 


Assuming  constant  L/D,  integrating,  and  using  a  reference  time,  tp,  to  provide 
the  constant  of  integration, 

^  "  •  Ycrf/DT  -  ^o) 

Defining  a  nondimens ional  time  variable  (Ref.  28), 

=  Vc(l/d)  -  t)  =  yff:  ^(l/d)  ^ 
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Eq.  (86)  can  be  written 

u  =  ^  =  tanh  T  or  U  =  Y^.  tanh  T  (83) 

*c 

The  constant  to  is  selected  to  fit  the  initial  conditions  pertinent  to  the 
problem  at  hand. 

The  dynamic  pressiire  may  be  found  by  noting  that  only  small  flight 
path  angles  are  permissible  if  L/d  tan  7  <C  1 ,  so  that 


cos  7  =  ^  (1  -  il^)  =  (1  -  u2)  4  1 


and 


1  ‘ 


^  (1  -  tanh^  T) 


=  sech^  T 

Cjp 


(89) 


The  density  ratio  p/p^  is 


^L^Ppfc 


Po  P„u2 


csch^  T 


-Ph 


and,  since  p  =  p  e  ^  ,  the  altitude  is 


(90) 


(91) 


The  flight  path  angle  has,  to  this  point,  been  specified  only  to 
the  extent  that  cos  7  must  approximate  unity  and  that  the  assumptions  of  Table  HI 
hold.  The  actual  angle  may  be  found  by  using  Eqs.  (87)  and  (91)  to  obtain  dh/dt 
and  then  noting  that  tan  7  =  h/U.  The  result  is 


1 

tan  7 


£l£  k 

2  D 


tanh^T 


L  ? 
"  2  D  ^ 


(92) 
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Again  y  is  always  negative  and  further 


Assigning  a  maximum  value  of  L/D  tan  y  -  O.O5,  for  example,  corresponds  approx¬ 
imately  (using  a  nominal  value  of  Pfc"  9^^  requiring  It  to  be  greater  than 
0.21  —  not  a  very  stringent  limitation. 

Plots  which  Illustrate  the  variation  of  U ,  h ,  and  q  with  the  non- 
dimensional  time  parameter,  T,  are  given  in  Figixre  25. 


Figure  23.  Time  Variation  of  U, h  ,  euid  q.  Gliding  Flight 


60 


5-  Skipping  Flight 

For  the  aBSun5>tlonB  pertinent  to  the  skip  situation,  the  simplified 
equations  of  motion  are 


U  =  - 


cos  y 


*  _  Pog~^*^U^  / 

2  cob2  r 


sin  7  “  2  cos  y 


Also,  by  differentiating  the  familiar  W  =  -  Utan  y, 


W  =  -Uy  sec2  7  -  U  tan  y 


Substituting  U  and  W  from  Eqs.  (95)  > 


-Bh  o  "Pf*  2  j 

2  Poe'^U^  /L  \  Pq®  U""  sin  7 

=  TTTT - i-  ( D  TT^ - i“ 

2  C062  7  \  >2  cob2  y 


(c^s) 


7  sec2  7  = 


(ctf) 


COB  7 


Assuming  constant  L/d,  integrating,  and  setting  the  initial  conditions  U^,  y^ 


^  =  2  (tan  7p  -  tan  y) 


iL  _  _(D/L)(tan  y^  -  tan  y) 


The  flight  path  angle  variation  with  altitude  may  be  found  by 
recognizing  that  h  =  U  tan  y.  Then 
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(99) 


Pq  Pq 


(sec 


sec  7) 


and  the  dynamic  pressure,  q,,  becomes 

.2 


_eu. 


2  cos^  7 


=  P 


/  m  \ (sec  7o  -  sec  7)  p  2(d/L) (tan  7n  -  tan  7) 
^cls)  ::z5“: 


"5" 

cos*^  7 


(100) 


The  variation  of  altitude,  density,  dynamic  pressure,  and  velocity 
as  a  function  of  7,  for  L/D  =  2,  is  shown  in  Figure  2h. 


o-i'ir  2 

ally  ^ 

r  ■ 

313  i 


X(deg) 

note  subscript  y  indicotes  condition  of  y  =0 


Figure  2k.  Variation  of  U,  h,  q,  and  p  for  Typical  Skipping  Flight 
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The  three  limiting  types  of  flight  considered  above  are  of  interest  in 
the  present  context  primarily  because  they  represent  siii5)le  particle  dynamics 
situations  for  which  the  equations  may  be  solved.  Table  IV  summarizes  the  im¬ 
portant  results. 

TABLE  IV 

SUMMARY  OF  DYNAMIC  FLIGHT  PARAMETERS  FOR 
READILY  SOLVABLE  CASES  OF  AEROSPACE  VEHICLE  PARTICLE  DYNAMICS 


SITUATION 

CHARACmiSnC  PARAKETER 

Skip 

Borltoatal  velocity, U 

Unh  T 

I _ 

y  ^(D/I.)(t.n  To  -  t«”  r) 

Flight  pftth  Angle,  7 

To 

Dyua&lc  preeiure,  q 

0  /  ,  V  t«n  ro  r  (f  . 

(cds)  CO.  ro  J 

Density,  p 

/_5_\  ^  ’’o  -t 

^  \cds)  CO.  To 

In  a  more  practical  vein,  the  three  limiting  cases  are  of  great  interest 
from  the  standpoint  of  vehicle  performance.  The  power-off  flight  performance 
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problem  may  be  considered  to  involve  fundamentally  the  conversion  of  initial 
velocity  into  range  while  maintaining  the  vehicle  structural  Integrity  in  the 
presence  of  critical  aerodynamic  heating  and  loading  effects.  This  problem, 
which  has  been  extensively  studied  and  documented  in  recent  years,  is  beyond  the 
scope  of  this  report  (see  Ref.  25-27). 


65 


CHAPTER  V 
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_  2GM 
dg  -  --^ir 


-2g 


dr 


dh 

r 


dh  =  -dW  =  tan  ydU  +  U  sec^  ydy 


where  the  single  F  is  the  ideal  power-off  constant  speed  glide  angle,  tan  F  =  -D/L, 
and  R  is  the  resultant  aerodynamic  force,  R  =  Vl^  +  D^* 

The  linearized  equations  of  perturbed  motion  then  become: 


lO*  [^t»n  7  ♦  (r  -r)]4u|  -  tm  7  ♦  e  I  *1"  (r  -  njdh  [a  .In  (7  -  D  tM  7  ♦  ».  (7  -D]  ♦ tec27|47  -  0 


[t  *  n  “*  -jdi; « ^  ♦  e ;  =0.  (7 


.  un  7<1U  *11*  *  0 


►  5  Ja  CO.  (7  -  D  iM  7  -  .In  (7  -r)]d7 


Eqs.  (102)  may  be  interpreted  to  represent  small  perturbed  motions  either  about 
flight  operating  points,  or  from  initial  flight  conditions.  They  are  precise  in 
the  differential  form  shown.  Because  of  the  nature  of  total  differentials,  the 
coefficients  enclosed  in  square  brackets  are,  at  worst,  functions  of  time  only, 
so  the  equations  are  linear.  Unfortunately,  unless  restrictions  are  applied  to 
keep  the  coefficients  constant,  solution  will  be  just  about  as  difficult  as  for 
the  original  unperturbed  set.  In  general,  then,  Eqs.  (102)  must  be  either 

a.  Linear,  but  with  time  varying  coefficients 

b.  Linear,  with  constant  coefficients,  but  restricted  as  to  the 
time  span  over  which  the  solution  is  valid 

For  the  equations  to  take  on  a  constant-coefficient  form  requires  either: 
the  practical  existence  of  a  condition  of  equilibrium,  where  the  operating  point 
conditions  are  described  by  constant  values  of  the  dependent  variables  U  and  W; 
or  the  assumption  that  the  coefficients  are  evaluated  at  some  starting  point 
(t  =  O) .  In  either  case  the  equations  will  be  good  approximations  to  the  motion 


(102) 
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only  as  long  as  the  perturbed  motions  remain  small  relative  to  the  total  values 
of  the  dependent  variables. 

Defining  a  steady- state  value  of  W  other  than  zero  implies,  unfortunately, 
that  there  is  a  steady-state  value  of  h,  since  W  =  -n.  The  perturbed  altitude, 
dh,  will  then  have  a  component,  -W^jt,  which  varies  with  time,  and  the  linearized 
equations  will  be  time  varying.  This  fundamental  problem  is  not  new.  It  was  en- * 
co\antered  in  aeronautical  stability  and  control  many  years  ago  (Ref.  29-32)  for 
diving  and  climbing  stability  calculations  in  which  density  variations  were  con¬ 
sidered.  In  those  Instances  a  steady-state  operating  point  corresponding  to 
7  =  0,  U  =  constant  resulted  in  tractable  and  meaningful  equations  of  motion.  In 
the  present  case  density  gradient  is  still,  together  with  gravitational  gradient 
effects  ^see  the  coefficients  of  dh  in  Eqs.  (102)3,  a  basic  source  of  difficulty. 
Now,  however,  the  glide  case  for  y  =  0  is  not  characterized  for  very  long  periods 
of  time  by  nearly  constant  values  of  U  (see  Fig\ire  23) •  Nevertheless,  the  only 
pertinent  operating  point  (or  initial  point)  conditions  yielding  solvable  equa¬ 
tions  are,  as  before,  those  for  y  =  0,  U  =  constant.  The  results  so  obtained  are 
still  of  considerable  physical  significance,  and  can  be  used  to  describe  the  par¬ 
ticle  motions  by  a  series  of  "frozen"  system  solutions.  After  setting  the  oper¬ 
ating  point  flight  path  angle,  y,  to  zero  and  expressing  L,  D  in  terms  of  basic 
aerodynamic  coefficients,  Eqs.  (102)  become: 


Eliminating  the  perturbed  flight  path  angle  gives. 
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If  the  operating  point  altitude  and  speed  are  assumed  to  he  constant,  the 
coefficients  confuted  with  operating  point  values  (denoted  hy  a  zero  suhscript) 
will  also  he  constant.  Identifying  the  differentials  with  small  perturbations , 
i.e.,  dU  =  u,  dh  =  h,  etc.,  and  Laplace  transforming  with  zero  initial  conditions, 
yields 


V  ' 

\  K 

mUoi  ^ 

+ 

_Uu^^  To- 

h  -  — _ 

=  0 


(105) 


/Uo 

lo  \ 

'  2  °° 

/ui 

2go 

3Lo\ 

k 

+ 

- 

f 

k' 

—  jj 

Although  Egs.  (105)  can  easily  he  studied  as  written,  it  is  more  pertinent, 
in  view  of  the  preceding  chapters,  to  specialize  them  into  extra-atmospheric  and 
atmospheric  cases.  Outside  the  atmosphere  (Lq  =  Dq  =  0)  Eqs.  (105)  reduce  to  the 
elementary  form 


(106) 


These  correspond  to  small  perturbations  about  essentially  circular  orbits, 
characteristic  equation  is 


A(s) 


The 


(107) 


The  last  result,  involving  the  mean  orbital  frequency  follows  directly  from 
Eqs.  (48)  and  (49).  Thus,  the  period  of  the  low  frequency  particle  motions  of  the 
vehicle  approaches  the  orbital  period  as  the  vehicle  leaves  the  atmosphere. 

B.  LINEARIZED  EQUATIONS  FOR  ATMOSPHERIC  FLIGHT 

For  flight  within  the  atmosphere  the  major  region  of  Interest  corresponds  to 
the  condition 
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For  such  situations  theUW/r  term  may  be  neglected,  as  shown  in  the  preceding 
chapter.  Also, 

,/u2\  2U  /  1  dr/r\  .  2U 

”2  dip/  -  -T-dU 

The  use  of  these  assumptions,  in  evaluating  the  total  differentials  connected 
with  Eq.  (101 ),  reduces  Eqs.  (l04)  to 

[t  *  3  (¥)]'!“}  -  *  3  (t)  “  <'“»> 

The  relative  magnitudes  of  the  terms  in  the  coefficient  of  dh  in  Eq.  (109)  are, 
making  use  of  Eq.  (82)  for  L/mg, 


is 

'  =  —  ^ 

pr(l  -ii^) 


Since  Pr  is  approximately  900^  the  ratio  is  small  with  respect  to  1  (about  5/S) 
for  values  of  u<  O.98.  Then  Eq.  (109)  becomes 


{  =  [7*3  1'*'^  *  3  (^)  O 


and  the  specialization  of  Eqs.  (105)  for  the  atmospheric  flight  situation  is 
given  by. 


(150) 


The  characteristic  equation  for  the  set  is,  neglecting  2  relative  to  Pfo, 

Before  examining  this  result,  a  similar  development,  assuming, 
however,  the  existence  of  a  constant  thrust,  Tq  =  Dq,  which  permits 
the  operating  point  to  he  in  complete  equilihrium,  will  he  tinder- 
taken.  Accordingly,  the  sets  of  equations,  Eqs.  (101)  and  (107), 
(neglecting  the  defining  equation  for  y)  are  modified  hy  the  addi¬ 
tion  of  the  following  Indicated  thrust  terms: 


UW 

u  -  ^  = 


To 

+  —  cos  y 


w+  ^  = 


-  —  sin  7 
m 


dU  + 


+  —  sin  7d7 
m  '  ' 


+  —  cos  7d7 
zn 


dU  + 


+  ^d7 

m 


dU  •+ 


(104,  108)' 


(104,  109)' 


+  ST 


(105,  106,  110)' 


Recognizing  that  Do/mh^ 
Eqs.  (no)  becomes 


=  To/mU^j,  the  complete  revised  set  of 
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and  the  characteristic  equation,  again  neglecting  2  relative  to  pr^,  is 

The  classical  two  degrees  of  freedom  siibsonic  phugoid  oscillation 
has  a  characteristic  equation  (e.g..  Ref.  3)^ 


=  0 

“^o  Ug  V  “  / 


derived  for  conditions  where  Tq  =  Dq,  lift  and  weight  are  equal,  i.e., 
Lo/m  =  6,  and  density  gradient  is  neglected,  i.e.,  p  =  0.  With  these 
same  restrictions,  Eq.  (113)  reduces  to  Eq.  (11^),  indj.cating  that  the 
first  three  terms  of  Eqs.  (113)  and  (111)  correspond  to  a  mode  akin  to 
the  familiar  phugoid. 

The  correspondence  noted  in  the  above  digression  can  be  further  clar¬ 
ified  by  approximately  factoring  the  cubic  equations,  (ill)  and  (113),  in 
llteial  terms,  i.e.. 


s3  +  As^  +  Bs  +  C  = 

(•  *  i)| 

s2  + 

.  (a  -  1)  .  b] 

which  is  valid  when 

AC  ^  1 
-5  <  1 

and 

£ 

b5 

1 

(11.^) 


For  both  equations,  the  B  coefficient  is  dominated,  for  u^  greater  than  about 
0.15,  by  the  term 

PUo  W  P«'o  ..  . 

2S  "  2  T  _  ^2  ^  ' 

so  that  3  =  ^  =  g^  ^  =  gP(l  -  u^)  (116) 


and,  neglecting  the  3:2  differential  between  the  A  coefficients  in  Eqs.  (ill) 
and  (11 3), 
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and 


Ib3 


/^Uq  Do\^ 

^  \  '■o  ^  4  u- 

,5  "  P^o(l)  1  -  ^  ^ 


ttS 


m 


With  the  conditions  of  validity  thus  shown  to  he  satisfied,  the  characteristic 
equations  have  the  approximate  factors 


A(s)  =  ^s  +  +  2^p“pB  + 


where 


2U„ 


(l/d) 
Lo 


2Vg/ro  _ 
L/D  '^o 


u|  =  P  ^  =  gP(l  -  ^o) 


2^p% 


2D. 


®Uq(1  -  u^) 


2g 

Uo(L/D)  "  Uo(L/d) 


for  To  =  Dq 


(117) 


^'’o  (d)^  V’  "  “o 

("D 


o 


for  To  =  Do 


for  To  =  0 


(5)^0 

The  variation  of  these  parameters  for  To  =  0  is  shown  in  Figure  25 


The  physical  forces  involved  in  each  of 
the  modes  can  he  studied  hy  constructing  time 
vector  diagrams  (Ref.  35  and  3^)  for  the  typi¬ 
cal  zero-thrust  example  given  helow: 


=  0-95 

2.0 


o 

L 

D 


To  =  21  .1  X  10'= 

P  =  23,500-'' 
g  =  32.2 


2Do 

mUb 


mU„ 


=  0.127  X  10' 


-5 


U„  = 


0.1270  X  10"5 
2li,750 


—  =  1 .175  X  10'^ 
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Figure  25.  Variation  of  Charac¬ 
teristic  Equation 
Approximate  Factors 


Eqs.  (no)  then  become 


(s  +  0.1270  X  10"^)u  +  (0.1270  X  lO'^s  -  0.067  X  10"5)h  =  0 
2.6o4  X  10“ -  (s^  +  0.0635  X  10“^s  +  0.1358  X  10"5)h  ==  0 


and  the  exact  factors  of  the  characteristic  equation  are 

(s  -  0.001 165) (s2  +  0.001356s  +  0.0001357) 
or,  (s  -  0-001 165) [s^  +  2(0.0586) (0.01 165)3  +  (0.01165)^] 

The  approximate  factors,  corresponding  to  Eqs.  (117),  have  the  values 

(s  -  0.001175)(s^  +  0.001363s  +  0.0001336) 
or,  (s  -  0.001175)[s2  +  2(0.0590)(0.01156)s  +  (0.01156)^ 


and  are  Seen  to  be  within  a  few  percent  of  the  exact  factors. 

The  time  vector  diagrams  of  the  lift  and  drag  equations  for  both  the  phugoid 
and  aperiodic  modes  are  shown  in  Figure  26.  Of  particular  interest  is  the  fact 
that  the  important  drag  equation  components  are  the  same  (i.e.,  h  and  li)  for 

both  modes.  Therefore  a  study  of  the  major  components  of  the  lift  equation  should 
reveal  the  essential  difference  in  the  physical  forces  active  in  each  mode.  The 
sin^jlified  lift  equations,  based  on  the  vector  polygons  of  Figure  26,  are 


Phugoid  mode; 
Aperiodic  mode: 


0 


It  appears,  therefore,  that  a  downward  change  in  height  is  accompanied  by  an 
upward  acceleration  giving  rise  to  the  phugoid,  and  by  a  decrease  in  forward 
speed  giving  rise  to  the  aperiodic  divergence.  For  the  phugoid,  all  the  es¬ 
sential  dynamics  (except  for  minor  damping  effects)  are  contained  in  the  sim¬ 
plified  lift  equation  which  directly  yields  the  approximate  phugoid  frequency 
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previously  derived  [Eq.  (117)1*  for  the  aperiodic  motion,  the  simplified  drag 
equation  must  also  be  used  to  obtain  the  complete  dynamic  picture,  the  slrplified 


set  then  being 


where  Lo/mLUj  has  been  neglected  with  respect  to  Uq/tq  consistent  with  the  con¬ 
ditions  chosen  for  this  example.  This  set  yields  the  time  constant,  Tp,  given 
in  Eqs.  (l IT) . 

On  the  basis  of  the  foregoing  exploration,  the  aperiodic  mode  contains 
most  of  the  elements  of  the  average  motion,  as  given,  for  example,  by  the  glid¬ 
ing  flight  equations  of  the  previous  article  and  pictured  in  Figure  23. 

The  phugoid  oscillation  is  associated  with  the  "boimcing"  commonly  encountered 
in  re-entry  maneuvers,  and  is  therefore  related  to  the  skipping  flight  situation 
treated  previously.  In  either  mode  the  major  aerodynamic  influences  act  only  by 
virtue  of  the  density  gradient  effect,  g.  This  explains  why  the  addition  of  a 
constant  thrust  has  so  small,  an  influence  on  the  characteristic  motions.  This 
would  not  be  the  case  if  thi’ust  were  assumed  to  vary  with  atmospheric  density. 


The  above  observations  must  be  restricted  to  time  spans  consistent  with  valid 
solutions  to  the  constsint  coefficient  equations  of  motion  utilized.  For  some  in¬ 
dication  as  to  how  long  such  times  may  be,  and  in  line  with  the  preceding  discus¬ 


sion,  assume  that  the  mean  time  history  of  the  vehicle-particle  is  given  by  the 
gliding  flight  solution.  Then  the  average  dynamic  pressure,  q,  will  be,  from 
Eq.  (89), 


=  ^  (1 
ClS 


Taking  the  time  derivative. 


and  utilizing  Eq.  (85)  for  du/dt. 
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Figure  26.  Time  Vector  Diagrams 


Since,  from  Eg.  (llT)>  the  approximate  value  of  o:^  is  proportional  to  g, 


1_  ^ 
2  g 


T/d 


udt 


The  time  span  for  which  the  constant  coefficient  or  "frozen"  solution  for 
will  hold  within  10^^  is  then 


At  -  0.10  Jiiii  i  52h/2l 

1  S  “  “o  “o 

Expressed  as  a  fraction  of  the  phugoid  period,  Pp  =  2rt/a^, 


^  =  Q-1Q  L  ~  ^  ^  0.48 


2it 


uo 


These  relationships,  plotted  in  Figure  27,  indicate  that  the  phugoid  freguency 


Figure  27-  Time  Span,  At,  for  "Valid"  Constant  Coefficient  Solutions 
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will  remain  within  1 0jo  of  its  "frozen"  value  for  times,  proportional  to  L/D, 
which  range  from  about  one-tenth  the  phugoid  period  at  Uq  =  O.98  (and  L/D  =  1 ) 
to  values  greater  than  the  period  for  u^  less  than  about  0.45.  While  a  solution 
good  for  only  one-tenth  the  i»eriod  of  the  resulting  motion  might  be  considered 
unreasonably  poor,  the  fact  that  the  actual  time  involved  is  of  the  order  of  1 00 
seconds  is  a  hi^ly  mitigating  circumstance.  For  example,  most  flight  control 
systems  will  invariably  have  response  times  at  least  an  order  of  magnitude  faster, 
and  analysis  of  such  systems  will  not  be  affected  by  a  lOO-second  time  limit. 

Also,  because  of  the  very  slow  time  variations  involved,  the  "frozen"  solutions 
can  be  used,  with  appropriate  initial  conditions,  as  point-by-point  approxima¬ 
tions  to  the  real  time-varying  situation.  Reference  28  presents  comparisons  of 
the  approximate  value  of  of  Eq.  (11?)  (differently  derived)  with  values  of ■  the 
skipping  frequency  obtained  from  complete  time-varying  numerical  computations. 

In  this  comparison  the  approximate  value  of  oUp  is  calculated  as  if  the  operating 
point  were  being  continuously  varied,  i.e.,  u  replaces  u^  in  Eqs.  (11 7)*  The 
agreement  is  excellent,  the  approximate  values  being  lower  than  those  determined 
numerically  by  less  than  about  5/^* 
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APPENDIX 


SUMMABY  OF  AXIS  SYSTEMS  FOR  AEROSPACE  VEHICLE 
STAETUTY  AJID  COOTROL 


This  appendix  summarizes  the  major  types  of  axis  systems  required  for  an 
adequate  analytical  description  of  aerospace  vehicles.  It  contains  a  series  of 
axis  system  definitions  and  illustrations,  a  simimary  table  of  these  reference 
systems  and  associated  nomenclature,  and  a  chart  which  summarizes  the  most  com¬ 
mon  transformations  between  systems. 

The  systems  presented  are  intended  to  be  adequate  only  for  aerospace  stabil¬ 
ity  and  control  situations  where  the  principal  gravitational  effects  upon  the 

* 

vehicle  rigid  body  dynamics  are  due  to  the  earth's  field.  Inertial  space  is  ac¬ 
cordingly  defined  as  a  geocentric  equatorial  system  referred  to  the  First  Point 
of  Arles.  In  this  situation,  the  effects  of  other  bodies  are  considered  as  forc¬ 
ing  functions.  For  regions  where  it  is  essential  that  forces  due  to  the  sun  and/ 
or  moon  be  considered  as  dependent  on  vehicle  attitude  (i.e.,  producing  terms 
which  would  appear  in  the  homogeneous  equations  of  motion),  additional  systems 
are  required.  For  the  sun  and  moon  these  would  be,  respectively,  a  heliocentric 
ecliptic  system  and  a  selenocentric  equatorial  system.  If  these  systems  were 
added  to  those  shown,  the  definition  of  inertial  space  would  be  shifted  to  the 
heliocentric  ecliptic  system. 

The  ground  rules  used  in  defining  the  systems  and  nomenclature  shown  are; 

1  .  The  systems  defined  must  be  sufficient  for  use  in  describing 

aerospace  vehicle  dynamic  situations  in  which  approximate  dynamics 
of  the  vehicle  plus  some  controller  are  of  major  interest.  For 
some  purposes,  such  as  early  boost  or  late  re-entry  phases,  the 
axes  involved  in  conventional  aerodynamic  stability  and  control 
(as  used  here,  XYZ,  Xp  Yp  Z^,  Xj^  Yj^  and  XYZ)  plus  a  geocentric 
equatorial  inertial  system  and,  perhaps,  attitude  sensor  reference 
axes,  wo\ild  be  sufficient.  In  other  phases,  however,  the  geocentric 
radius  is  not  necessarily  an  adequate  vertical  reference  (Ref.  35). 


In  a  more  general  sense,  the  systems  are  suitable  for  situations  where  a  solitary 
principal  gravitating  mass  is  present. 
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This  point  is  emphasized  in  Figure  A-1 ,  which  illustrates  some  ver¬ 
tical  reference  possibilities.  Directional  references  are  similarly 


GEOCENTRIC: 

n.  INSTANTANEOUS:  Gnocentric  rndiun  throucli  the 
vnhicln:  line  .  (^  auTm) 

b.  MEAN  NOTION:  Gnocentric  rndiun  throufh  the  point 

which  Bovns  with  thn  vehicle  aenn 
Bot  i  on . 

c.  CONSTANT  VELOCITY;  Geocentric  rndiuo  throufh  the 

point  which  would  Bove  with 

the  vwhicln  Been  Botion  if 

there  were  no  orbltnl  rnirnssion. 


gravitational: 

Linn  through  the  vehicle  norael  to  the  nsoothed 
ellipnoldnl  nppr oe iBe 1 1  on  to  the  plnnet  grn v i te 1 1 onn 1 
nauipotcntlnl  nurfnce  passing  through  the  vehicle; 
line  GV. 

GEOGRAPHICAL: 

Line  through  the  vehlclH  noraal  to  the  1 nte rna t  1  one  1 
Ellipsoid;  line  KV. 

HORIZON  SCANNER: 

Linn  through  the  vehicle  through  the  centroid  of  surfnce 
bounded  by  the  horizon;  line  Hv. 


Figure  A-1  .  Vertical  Reference  Possibilities 


nonunique.  Accordingly,  several  additional  systems  are  required 
to  provide  for  adternate  attitude  references  and/or  departure  of 
the  vehicle  from  a  Keplerian  trajectory. 

2.  The  notation  advanced  should.  Insofar  as  possible,  be  compatible 
with  both  astronautical  and  aeronautical  usage . 

5.  Both  the  systems  chosen  and  the  notation  used  should  reduce  to 

those  used  in  aeronautical  stability  and  control  as  one  limit.  As 
can  be  detected  by  a  detailed  examination  of  the  figures,  this  cri¬ 
terion  makes  necessary  an  awkward-appearing  choice  in  the  definition 
of  the  geocentric  equatorial  inertial  system,  i.e.,  the  axis  is 
directed  toward  the  First  Point  of  Aries.  This  condition  derives 
from  the  decision  that  the  vehicle  body  axis  system  should  corre¬ 
spond  to  that  conventionally  used  on  aerodynamic  vehicles,  viz.,  X 
"forward,"  Y  out  the  right  wing,  and  Z  "down."  Starting  with  this 
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latter  system,  and  a  west  to  east  orbit,  and  working  backward 
makes  the  result  noted  inevitable.  The  reason,  of  course,  for 
choosing  the  body  axis  system  as  a  starting  point  is  to  preserve 
unchanged  the  extensive  notation  and  physical  appreciation  already 
built  up  for  aerodynamic  vehicles. 

4-.  When  possible,  the  rotatioMil  sequences  involved  in  going  from 
one  system  to  another  should  be  similar.  Thus,  the  various  axes 
in  the  vehicle  are  derived  by  similar  rotations  from  the  body- 
fixed  reference  axes. 

5.  The  number  of  symbols  used  should  be  minimized,  consistent  with 
clarity.  Thus,  the  unit  vectors  i  j.  k>  with  appropriate  sub¬ 
scripts,  are  losed  for  all  vehicle -centered  systems  which  do  not 
contain  some  externally  based  vertical  reference  possibility  as 
an  axis.  Similarly,  the  Euler  angles  relating  all  vehicle-cen¬ 
tered  systems  are  taken  as  ♦,  8,  and  ♦,  with  appropriate  sub¬ 
scripts.  This  selection  preserves  well  imderstood  physical 
notions  of  yaw,  pitch,  and  roll  motions  of  one  system  relative 
to  another. 

6.  The  symbols  used  shoxild  be  easily  reproducible.  This  criterion 
has  been  met  as  long  as  one  has  handy  a  typewriter,  standard 
Leroy  gothic  alphabet,  and  someone  capable  of  writing  X,  Y,  and 
Z  in  script. 

Of  all  these  requirements  and  desires,  the  third  is  unquestionably  the  most 
important  because  aerospace  vehicles  are  also  aircraft.  The  axis  systems  se¬ 
lected  to  satisfy  the  third  point,  however,  need  not  be  confined  to  those  defined 
here .  M.  Abzug,  also  recognizing  the  logic  of  compatibility,  has  thoroughly  exam¬ 
ined  most  of  the  possibilities  (Ref.  36).  He  notes  that  some  awkward  features, 
such  as  the  direction  of  the axis  toward  T,  negative  rotation,  etc.  are  inev¬ 
itable.  Some  degree  of  arbitrariness  also  exists,  such  as  easterly  or  westerly 
headings,  etc.  Therefore  the  axes  chosen  here  are  not  all  uniquely  defined  by 
the  requirements  and  desires  (except  for  the  body-centered  systems) .  However, 
the  end  results  presented  here  exe  reasonable  and  should  provide  at  least  a  basis 
for  constructive  criticism. 
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The  various  systems,  transformations,  and  summary  tables  follow.  Figure 
A-2  shows  the  basic  steps  in  the  development  sequence  from  inertied  space  to  the 
fundamental  body-centered  set,  XYE,  via  the  trajectory  defining  quantities  Afl, 

L,  w,  and  V.  This  entire  sequence  can  be  specialized  to  a  minimum  of  two  sys¬ 
tems,  as  shown  in  Figure  A-3,  if  oblateness  is  ignored  or  negligible.  An  alter¬ 
nate  route  from  inertial  space  to  3CY2  is  illustrated  in  Figure  A-4.  This  se¬ 
quence  introduces  latitude,  longitude,  and  northerly  heading  angle  in  lieu  of  the 
trajectory  elements. 

The  systems  noted  above  are  the  only  ones  actually  used  in  this  report. 

When  rigid  body  dynamics  are  considered,  however,  further  systems  are  needed. 

The  fundamental,  body- fixed,  set  XYZ  is  shown  in  Figure  A-5.  Its  relationship 
with  the  body- centered  XYZ  system  is  shown  in  Figure  A-6.  The  rotational  se¬ 
quences  between  XY2  and  XYZ  are  summarized  in  Table  A-I.  Similar  sequences  can  be 
used  to  relate  principal  axes,  Xp  Yp  Zp,  moving  element  axes  Xj^  Y^  auxiliary 
attitude  axes,  xyz,  etc.  as  required  with  either  XYZ  or  XYZ. 

Finally,  Table  A-II  recapitulates  the  various  sets  and  relationships  and 
indicates  a  tentative  nomenclature  for  unit  vectors  and  resolution  angles. 


Figure  A-5-  Simplified  Situation  for  Homogeneous 
Spherical  Central  Body 
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NORTH  POLE 


♦ 


(Ascending  Node  of  Sun) 


Ascending 
Node  of  Vehicle 


ARIES  DIRECTED 
GEOCENTRIC  AXES 


ASCENDING  -  NODE  -  DIRECTED 
GEOCENTRIC  AXES 


Derived  from  rotation 

about  %  through  ,  the  longi¬ 
tude  of  the  vehicle  trajectory 
ascending  node. 


Figure  A-2.  Basic  Axis  System  Development 
(Homogeneous  Spherical  Earth) 
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Perigee 


Trajectory 

Reference 

Plane 


^Equatorial  Plane 

VEHICLE- CENTERED 
)CENTRICALLY-DIRECTED 
AXES  (XYZ) 

Lved  by  rotation 

through  V,  the  true  anom- 
;  and  then  translation  along 
to  the  vehicle  center-of-mass . 


PHYSICAL  SITUATION 

Aries -Directed  Geocentric  System 

XYE*  Northerly-  and  Qeocentrlcally-Dlreoted 
Vehicle-Centered  System 

X'  In  Direction  of  Increasing  Latitude 

Y'  In  Direction  of  Increasing  Longitxide 

Z'  Directed  Down  Geocentric  Radius 

XYE  Vehicle-Centered  Geocentrically-Directed 
System 

Figure  A-4-.  Supplementary  Relationships  Between  Vehicle-Centered 
Geocentrically  Directed  and  Aries -Directed  Geocentric 
Systems 


INTERMEDIATE  SYSTEM 

Derived  from  rotation  of 
about  %  through  the  longitude, 
A,  measured  from  the-^  axis. 


NORTHERLY-  AND  GEOCENTRICALLY - 
DIRECTED  VEHICLE -CENTERED 
SYSTEM  (X'rZ') 

Derived  from  hy  rotation 

about  Y"  through  the  latitude  angle 
X;  and  then  translation  along  -  Z*  to 
the  vehicle  center-of-mass. 


VEHICLE- CENTERED  GEOCENTRICALLY - 
DIRECTED  SYSTEM  (XYZ) 

Derived  from  X*Y^E*  by  rotation  about  E* 
throu^  the  northerly  based  heading  angle 

'I'n. 

Figure  k-h,  (Continued) 
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I  ACCELERATION 


ARROWS  INDICATE 
POSITIVE  SENSE 


^  ^  acceleration 
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T 
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R 


APPLIED 
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M 

N 


Figure  A-5.  Vehicle-Fixed  Axis  System  and  Notation 
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RELATIONSHIPS  BETWEEN  AXES  ORIENTED  TO  EACH  OTHER  VIA  AND<I> 
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Arxes-Dlrected  Geocentric  Syntei:i 


Ascendinp-Node-Di  re  etc  <1 


Geocentric  Axes 


Trajectory  Reference  Plane 
Geocentric  Axes 


Vehicle-Centered  Gt-oceutri  -aliy- 
■Oirected  System 


Northerly-  and  GcocentricaLly- 
Directed  Vehicle-Cv-nt'-r'-d  i'.yncei:; 


Vehicle-Fixed  R-feren'''-  AX'-s 


Vehicle  Principsl  Axi 


z'y'^ 


■%  y  y 


X  Y  Z 


V  Y  '■ 
'•p^iyp 


Attitude  Hel'crr'iiC' 


Moving  Element  A>;';n 
(Subscript  "ri"  ld(,'ntiri< 

inoviriG  <.-loi:ient) 


directed  along  north  pjolar  azio;-^( 
(Sun's  ascending  node  as  seen  from  th 


^^^^^rotated  about  through 
for  inertial  space  when  regression  of 


y ^  rotated  about  ^  through  t 
the  argument  of  perigee,  W  .  Approxi 
and  variation  of  the  line  of  apsides 


^‘p  rotated  about  through 
of -mass  (or  other  reference  point)  in 
stability  and  control  problems 


y  rotated  about  through  the 
rotated  about  through  the  vehicl 
(or  other  reference  point)  in  the  vel 
Occasionally  used  as  inertial  space  ; 

XYE  rotated  about  Z  through  yaw  ani 
X"Y"'i"  rotated  about  Y"  through  pitcl 
X'Y'Z'  rotated  about  X’  through  roll 

XYZ  rotated  about  Z  through  an  angle 
XpYpZp  rotated  about  Yp  through  an  a 
XpYpZ^  rotated  about  Xp  through  an  a 
Can  also  derive  from  XYE  by  rotatio 


XYZ  rotated  about  Z  through  an  angl 
xiyizi  rotated  about  yi  througli  an  ai 
^2X2^2  rotated  about  Xg  through  an  ai 
xyz  ordinarily  reduces  to  XYZ  when  1 

XY7,  rotated  about  7,  through  an  angle 
X|'jYj'|Z'(  rotated  about  Y"  through  an  a: 
X|'jY]'j!.^  rotated  about  X,'j  through  an  a: 


TABLE  A -II 

SUMMARY  OP  MAJOR  AXIS  OR  REFERENCE  SYSTEMS 


iS 

DESCRIPTION  and/or  DERIVATION 

ASSOCIATED  PHi 

1 

.^directed  along  north  polar  axis; directed  toT,  the  First  Point  of  Aries 
(Sun's  ascending  node  as  scon  from  the  earth) 

Basic  inertial  referen( 
vehicles 

/ 

^^^^^rotated  about  through  the  longitude  of  the  vehicle  ascending  node.  Approximation 

for  inertial  space  when  regression  of  the  line  of  nodes  is  negligible 

Rotation  of  line  of  no( 

2y 

^  ^ ^  rotated  about  ^  through  the  trajectory  coinclination,  it/2  -  i ;  then  about  "through 

the  argument  of  perigee,  W  •  Approximation  for  inertial  space  when  regression  of  the  line  of  nodes 

and  variation  of  the  line  of  apsides  arc  negligible 

Rotation  of  line  of  aps 

^ X  rotated  about  through  the  true  anomaly,  v;  origin  is  then  translated  to  the  center- 

of-mass  (or  other  reference  point)  in  the  vehicle.  Used  as  inertial  sx^ace  for  most  aerodynamic 
stability  and  control  problems 

Center  of  mass  movemeni 

1 

^  rotated  about  through  the  vehicle  longitude.  A;  (intermediate  system  X  Y  Z  )  then 

rotated  about  through  the  vehicle  latitude,  X;  origin  is  tlien  translated  to  the  center-of -mass 

(or  other  reference  point)  in  the  vehicle.  Related  to  XYE  via  the  northerly  heading  angle  SP'jj. 
Occasionally  used  as  inertial  space  in  aerodynamic  vehicle  automatic  control 

Implies  heading  refcrei 
torques  due  to  earth's 
longitude  variations  ii 

XY  Z  rotated  about  Z  through  yaw  angle  to  form  intermediate  .system  X"Y" 

X"Y"Z"  rotated  about  Y"  through  pitch  angle  0  to  form  intermediate  system  X'Y'  f 

X'Y'Z'  rotated  about  X'  through  roll  angle  to  form  XYZ 

Rigid  i'ody  laotions 

XYZ  rotated  about  Z  through  a,n  angle  ''1'’^  to  fonii  intermediate  system  XpYplp 

XpYpZp  rotated  about  Yp  through  an  angle  @p  to  form  intermediate  syssem  XpYpi.^,  i 

X^YpZ^  rotated  about  X^  through  an  angle  to  form  XpYpZp 

Can  also  derive  from  XYZ  by  rotation  through  similar  sequence 

XY  Z  rotated  about  Z  through  an  ang,!-'  to  form  tlic  intermediate  system  XjyjZ) 

xiyizi  rotated  about  yi  through  an  angle  to  form  tlie  iiitenncdiatc  systen. 

^2X2^2  rotated  about  xg  tlirough  an  angle  to  form  xyz 

xyz  ordinarily  reduces  to  XYZ  when  oblptniie-ss  and  attitude  scnsoi'  alignment  effects  are  removed 

Supplementary  system  ir 
for:  attitude  referenc 
XYZ.  More  than  one  £ 
quired,  e.g.,  xyz  coulc 
mean  motion,  and  nnothe 
define  attitude  sensor 

XYZ  rotated  about  Z  through  an  angle  to  form  the  intermediate  system  X^Yi”,:’./, 

^n^r?n  rotated  about  Y"  through  an  angle  0  -o  fonn  Llie  intcmiicdiatc  system  Xj',Y, 

X_)|Y]^Zj!j  rotated  about  X/j  through  an  angle  <l>[,  to  form  XnYuZn 

Supplementary  axes  reqi 
moving  parts  within  the 

TABLE  A -II 


F  MAJOR  AXIS  OR  REFERENCE  SYSTEMS 


DESCRIPTION  and/or  DERIVATION 

ASSOCIATED  PlfYSICAL  EFFECTS 

ar  axis; -^directed  toT,  the  First  Point  of  Aries 
een  from  the  earth) 

Basic  inertial  reference  frame  for  aerospace 
vehicles 

trough  the  longitude  of  the  vehicle  ascending  node  ,  Approxiination 

rgression  of  the  line  of  nodes  is  negligible 

Rotation  of  line  o’  nodes 

g  through  the  trajectory  coinclinatloii,  n/2  -  l;  then  about  yrp  Llirough 

V  .  Approximation  for  inertial  space  when  regression  of  the  line  of  nodes 

of  apsides  are  negligible 

Rotation  of  line  o;.'  apsides 

through  the  true  anomaly,  v;  origin  is  then  translated  to  the  crenter- 
ce  point)  in  the  vehicle.  Used  as  inertial  space  for  most  aerodytiairiie 

olems 

Center  of  mass  movemaent  in  trajectory  plane 

arough  the  vehicle  longitude,  A;  ( Intenuediato  system  X”Y”E^’)  th(  n 
the  vehicle  latitude,  X;  origin  is  then  translated  to  tlio  center-of-mass 
)  in  the  vehicle.  Related  to  XYE  via  tiie  northerly  headir-g  angle 
tial  space  in  aerodynamic  vehicle  autoi.'iatic  control 

Implii-'S  lieading  rei’f.-rence  (north);  gravity 
torques  due  to  eart)i'G  oblateness;  latitude - 
lC;:glt\jde  variatioi  s  ir.  gravitational  field 

ough  yaw  angle  'I'  to  form  intermediate  system  X'Y":.'," 
hrough  pitch  angle  0  to  form  iiaennedia  te  system  X'Y'  ' 

■hrough  roll  angle  to  foi'in  a'Y;-'. 

iiigid  :  ody  rjotions 

Lgh  an  angle  'I'p  to  form  interrae-dial  ■  sy.stem  X,"Yp '.p 
through  an  angle  0p  to  for:;,  inleritiediate-  syslei-  XpY|',.'.p 
;hrough  ar.  angle  4>p  to  form  X^YpZp 

1  by  rotation  through  similar  socpieiice 

rough  an  angle  to  form  the  intemicnlia If  systc;:;  X|y|:'.| 

through  an  angle  0,^  to  form  the  intennodiate  syster.  x-,y., 
through  an  angle  to  fon:i  xyz 

D  XYZ  when  oMa.'.eness  and  attitud--  sensor  alignment  effeern  are  r' :-.ov' ■  i 

;tippi<  mi'iitary  systei '  introduced  to  account 
for:  attitude  rcfel'cnce  axes  different  from 

XYZ.  Mor"  than  oii'-  Gin;h  set  is  often  re- 
quir'  d,  xys  ''C’  il'i  move  with  tlie  vehicle 

motion,  and  anooner  sot,  x'y's'  would 
defiin'  attitude  sensor  axes 

ugh  an  angl-'  to  \,\,(;  interned  iatr-  oyst(.-m  X.'JYJ', 

through  at£  angle  0  ..o  fonn  tiie  intcnned  late  system  y.'Y'.'' 

through  an  angle  to  fon:,  XriY,,Z,, 

riupplcmeiilai'y  axes  required  to  account  for 
movinc;  pai'ts  within  Idle  vcliicle 
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